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ANALYSIS OF PRESSURE-DROP FUNCTION IN RANKINE SPACE POWER 
BOILERS WITH DISCUSSION OF FLOW MALDISTRIBUTION IMPLICATIONS 
by Andrew A. Schoenberg and Donald R. Packe 
Lewis Research Center 

SUMMARY 

A problem in the development of single-pass boilers for Rankine space power sys- 
tems is the potential flow maldistribution of the working fluid among the parallel tubes of 
the boiler tube bundle. Flow maldistribution can result in undesirable liquid carryover 
to the turbine. This maldistribution can occur if the pressure drop of an individual tube 
is the same for more than one value of flow. This multivalued characteristic has long 
been recognized for boilers with low exit qualities; however, recent tests and the analysis 
presented herein show that a similar multivalued pressure-drop characteristic can also 
exist in boilers that operate with superheated vapor at the exit. 

The analysis was performed for counterflow liquid-heated single-pass boilers of the 
type used in Rankine space power systems such as SNAP-8. The analysis takes into 
account the nonuniform axial distribution of heat flux, variable geometry factors, and the 
variation of the lengths of the preheat, boiling, and superheat regions. Comparison with 
data from a single -tube water boiler indicates that the analytical model predicts the pres- 
sure drop and its dependence on various boiler -input variables with an average error of 
less than 10 percent. 

With the use of results generated from the analytical model, effects on the pressure- 
drop function caused by variations in some of the geometric and thermodynamic param- 
eters of the boiler were studied. From these studies, the following causes of the multi- 
valued characteristic were established: 

(1) The magnitude of the multivalued effect is strongly influenced by a boiler-tube 
insert, but the characteristic could exist even with an open tube. 

(2) The multivalued characteristic can exist for either a choked nozzle or a constant 
pressure at boiler discharge. 

Methods for eliminating the multivalued characteristic and thus the potential for flow 
maldistribution are evaluated. 


INTRODUCTION 


Rankine space power systems such as SNAP-8 (ref. 1) use a single-pass counterflow 
multitube boiler to vaporize the working fluid that drives the turboalternator. A problem 
with such boilers can be uneven division of the working -fluid flow among the parallel 
tubes. A particularly serious form of maldistribution, with attendant liquid carryover to 
the turbine, can result if pressure drop as a function of flow of the individual tubes is a 
multivalued function (i. e. , if the tube has the same pressure drop for more than one flow). 
Such multivalued pressure-drop functions, which are characterized by a negative slope of 
pressure drop as a function of flow for a range of working -fluid flow, have been studied 
for boiler tubes with low exit quality (refs. 2 to 6). However, for boiler tubes with an 
extensive superheat region (quality of 1. 0 at exit), the problem of multivalued pressure 
drop as a function of flow has received little attention. A recent test of a 19-tube counter- 
flow boiler (ref. 7) indicates that the pressure-drop function of these boilers can be 
multivalued even if superheated vapor is being generated. The nature of this character- 
istic has not been well understood, and indications are that a number of factors related 
to geometry, heat-flux distribution, and boundary conditions may be influential. 

The study reported herein, therefore, was conducted with the purpose of deriving a 
pressure-drop-against-flow function for a typical counterflow boiler tube in terms of 
dimensionless variables and parameters. This function would be useful for analyzing the 
factors that make the function multivalued. A further objective of this analysis was to 
determine effective methods of making the pressure-drop function single valued and, 
hence, of eliminating this potential for flow maldistribution in a multitube boiler. 

The analysis is limited to a counterflow liquid-heated tube-in-tube boiler and takes 
into account the axially nonuniform heat-flux distribution, variable -geometry effects, 
and the lengths of the preheat, boiling, and superheat regions. A check on the accuracy 
of the analysis was obtained by comparison with data of a single-tube water boiler. With 
the use of results generated from the analytical model, effects on the pressure-drop func- 
tion caused by some of the geometric and thermodynamic parameters of the boiler were 
examined, and from these effects the causes of the multivalued characteristic were 
established. Finally, methods for eliminating the multivalued characteristic and thus the 
potential for flow maldistribution were evaluated. 


DERIVATION OF BOILER EQUATIONS 

A typical tube in a single-pass boiler of a Rankine space power system is shown in 
figure 1. The heating fluid flows on the outside of the tube, counterflow to the working 
fluid on the inside of the tube. The wire-wound plug and the spiral wire shown are flow 
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swirl devices intended to promote heat transfer and phase separation by breaking up the 
larger liquid droplets as well as by centrifuging the denser liquid to the wall. 

For purposes of analysis, the boiler tube was divided into three regions: 

(1) The all-liquid region, where the working fluid is preheated to the saturation tem- 

perature 

(2) The boiling or two -phase region, where the liquid is vaporized 

(3) The gaseous region, where the vapor is superheated 

In all three regions only the contribution of friction to pressure drop was considered. 
Momentum pressure drop was neglected since it is negligible compared with friction 
pressure drop for the cases studied. For simplicity, the pressure gradient in the two- 
phase region was based on homogeneous flow, and a constant friction factor was deter- 
mined for a vapor quality of 1. The formulation of a pressure gradient included the vari- 
ation of velocity with local quality and took into account the nonuniform and variable heat- 
flux distribution and boiling length. The heat-transfer analysis was based on standard 
counterflow heat -exchanger equations. Constant, average heat-transfer coefficients were 
assumed in each of the three regions. For purposes of integration of the heat-flux equa- 
tion in the boiling region, a constant saturation temperature was also assumed. The 
validity of these assumptions is discussed during the development of the applicable equa- 
tions. 

In the first part of this section, the basic pressure-drop equations used in each 
region of the boiler are presented, and their dependence on the boiler heat flux and geom- 
etry is indicated. In the second part, the equations for the heat flux and the length of the 
preheat, boiling, and superheat regions are derived. Finally, the equations are synthe- 
sized, and the boiler -pressure-drop function is presented in dimensionless form. 

Basic Pressure-Drop Functions for the Three Boiler Regions 

For the single-phase (preheat and superheat) regions, a form of the standard tur- 
bulent friction-loss function (eq. (1)) was used: 

dP = fP fvV\- ( fp\ W » (1) 

di 4A C \2g j y gA 3j p 

(All symbols are defined in appendix A. ) The friction factor was assumed to incorporate 
the effect of the spiral wire and the effect of the bends in the spiral flow passages of the 
plug. The derivative dP/dZ was taken with respect to the axial distance along the tube 
so that in the plug region a multiplication factor dZ pl /d l must be used to account for the 
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extra length of the spiral path. An example of how these multipliers are evaluated is 
given in appendix B. K dP/d l is assumed to be constant in the single-phase regions, 
equation (1) can be integrated over the preheat and superheat lengths to give equa- 
tions (2) and (3), respectively: 


ap p 


f pgpV p 

Pw, l 


( 2 ) 


AP 


s 


f s g s w w*s 

\v 


(3) 


where 


g 


s 


Ps 

8gA* s 


g 


P 


p p i fegA 
8gA l ' 


The pressure drop in the two-phase region requires a more elaborate analysis, 
which includes not only the effect of the variation of length and geometry, but also the 
local variation in velocity resulting from quality changes. A correlation based on that 
proposed by W. L. Owens (ref. 8), and which assumes homogeneous flow, was used in 
this analysis. (A similar analysis for the case of uniform heat flux is given in refer- 
ence 9. ) In the analysis presented herein, several simplifications from Owens' analysis 
were made. The pressure drop due to change in elevation was neglected since space 
power boilers will generally operate in zero gravity. As mentioned earlier in this sec- 
tion, the momentum pressure drop was also neglected. Only the frictional component, 
of the following form, remains: 


dP 

dl 


jb_N' 

2Dg \A C 


_*_+(!. X )_L_ 

_ P W, V Pw, l_ 


(4) 
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The total pressure drop in the boiling region is expressed by 


AP b 



f b Sb w i 


W, V 


x + (l - x) 

- P w ,1 J 


dy 


(5) 


where 


g b =g p Osys(i pI -J p ) 

'b = S' > <*pl - 'p> 


( 6 ) 


Equation (5) is further simplified because the ratio of the densities p w / p w ^ is gen- 
erally small for the operating pressures of space power system boilers. Hence, it was 
assumed that a good approximation to the integral can be obtained by neglecting the term 
multiplied by p w y /p w i • The resultant integral takes on the simple form 

^ dy (I 

^W, V 

To permit closed-form integration of equation (6), it was further assumed that and 
P w y were constant so that they could be taken outside the integral, as shown in equa- 
tion (8): 



AP b = 


w. 


w 


W, V 



( 8 ) 


This simplified form can now be used to calculate the variation in pressure drop due to 
boiling length l b , geometry g b , and the local heat flux or quality x, as well as vapor 
density p w and total flow rate w w . The variation of x as a function of y remains 
to be determined from the heat-flux analysis. This analysis together with the determin- 
ation of the lengths of the three regions is presented in the next section. 
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Derivation of Heat Flux and Lengths of the Three Regions 


The heat-flux distribution in a counterflow boiler differs considerably from the more 
uniform distributions obtained with electrical or constant-temperature heat sources. 
Figure 2 shows a typical temperature and pressure profile in a counterflow boiler. The 
working fluid enters subcooled at a temperature T . and is rapidly heated to the satur- 
ation temperature T gat corresponding to the local saturation pressure P ga ^. In a 
counterflow boiler, the temperature difference T hb “ T sat» termed the "pinch-point 
temperature difference, " represents an important limitation on the amount of enthalpy 
available to vaporize the working fluid. The heating fluid, as shown by the upper curve, 
experiences little temperature change in the preheat region because of the small amount 
of enthalpy change required for preheating as compared to boiling. By far the greatest 
amount of heat-transfer and heating -fluid temperature drop occurs in the boiling region. 
This region of length begins inside the plug and normally extends over a large part 
of the tube length. The working -fluid temperature in this region remains at the satura- 
tion temperature, which decreases slightly because of the pressure drop in the boiling 
region. 

Of the total tube length, the portion L remaining after the preheat and boil- 
ing regions serves to superheat the vapor stream. For an actual boiler, the transition 
from boiling to superheating is not as distinct as indicated in figure 2 becasue of the non- 
equilibrium coexistnece of liquid droplets and superheated vapor over a considerable 
part of the tube length. For the purposes of the analysis, however, an abrupt transition 
was assumed. As in the preheat region, relatively little heat transfer occurs in the 
superheat region, as indicated by the small drop in heating -fluid temperature. The deri- 
vation of the equations for the three regions follows. 

The heat-transfer and fluid-exit temperatures in the preheating region can be deter- 
mined from the standard effectiveness relation for a counterflow heat exchanger (ref. 10): 

£ _ T sat- T wi _ l-exp[-N p (l-a p )L p ] 

P T hb- T wi l-a pe xp[.N p (l-a p ) L p ] 


where 






Note that N is a positive dimensionless parameter commonly referred to as the number 
of transfer units (NTU) of a heat exchanger which is herein based on the total tube length 


L t and the coefficient of heat transfer in the preheat region U 


tion (9) are known or will be specified except the preheat length L 
simplified as follows in order to solve for L 


All the terms in equa- 
Equation (9) can be 


_p. Given that Np is positive and that for 
most high-quality boilers the parameter a is much less than 1, the right side of equa- 
tion (9) can be expanded in a series: 


-N (1 - a )L " 

P P P_ ^ ^ _ 


1 - exp 
L-ap exp -N p (l-a p )L p 


(l-a p )exp[-Np(l-a p )Lp] -(l-a p )a p exp[-2N p (l -a p )L p ] . . . 

( 10 ) 


Neglecting all second order or higher terms in equation (10), substituting into equation (9), 
and solving for L result in 

Ir 


L 


P 


1 


- V N p 



(ii) 


This equation completes the derivation of the necessary equations for the preheat 
region. Next, the equations of the boiling region are derived. 

The heat flux in forced-convection boiling is a complex function of the flow regimes, 
temperatures, and surface conditions. A general discussion of these regimes can be 
found in references 2 and 4. Instead of attempting to predict the complex variations in 
boiling-side heat-transfer coefficient, an average constant overall coefficient was 
assumed for the calculation of heat flux in the boiling region. This assumption is justi- 
fied where the heating-fluid and wall -heat -transfer resistances are controlling. This 
case occurs in most space power boilers in the regions where contact boiling and, hence, 
high working -fluid heat-transfer coefficients are maintained (ref. 5). The region of con- 
tact boiling can be maintained to local qualities as high as 80 percent, particularly in 
boilers where the spiral wire helps to maintain annular flow (ref. 11). Thus, it is only 
in the last part of the boiling region that the boiling coefficient becomes small and, hence, 
controlling. Figure 3 illustrates the difference between the typical experimental overall 
heat-transfer coefficient and the coefficient assumed for the analysis. It also shows the 
resultant difference between typical experimental and computed temperature profiles of 
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the heating and working fluids. The saturation temperature was assumed to be constant 
at the value where the boiling region begins. 

With the assumptions of constant and uniform T g over the boiler length, the 
effectiveness relation (ref. 10) can be expressed as 


e b = 


T hs ~ T hb 
T hs " ^sat 


= 1 - exp(-N b L 1( ) 


( 12 ) 


where 


N b = 


V D o*t 

w h c h 


and 



The following conservation of energy for the whole boiling region can be applied to solve 
equation (12) for the boiling length in terms of specified variables: 

'Vh< T h S - T hb> = Wo < 13 > 


The quality is equal to 1 except in cases where the boiling length would exceed the avail- 
able tube length. This condition was included to account for the pressure drop in the 
boiler when the quality at exit is less than 1. Substituting the expression of equation (13) 
for the terms T hg - T hb in equation (12) and rearranging to solve for the boiling length 
result in 



(14) 


where 


€ b " 


w HX 
w w o 


w h c h< T hs - T sat> 
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The length of the superheat region L which is required to compute the pressure 

s 

drop in that region, is the length remaining after the preheat and boiling lengths, as 
given in dimensionless form in equation (15): 


L = 1 - (I. + L ) 


<Lb + L pXl'O 


L s = ° h, + s * l - 0 


(15) 


Note that + Lp cannot become greater than 1 (i. e. , exceed the length of the available 
tube). Hence, if the calculation of from equation (14) exceeds the length constraint, 
it means that the quality at the boiler tube exit is less than 1. Consequently, the outlet 
quality X Q which enters into the effectiveness expression must be adjusted so that the 
end of the boiling region corresponds to the end of the tube. With wet vapor at the tube 
exit, the appropriate temperature expressions are also adjusted, (i.e. , T hg - T h . and 
T wo = T ga j.). For normal boiler operation where superheat is present, the temperature 
expressions are given in the following equations: 


T hb = T ’ 


hs 


w w H w X o 

w h c h 


(16) 


T hs = 


hi 


WO 


sat' 


(17) 


where 


w c 

= W w,v 

s " Vh 


wo 


= T. 


hi 




(18) 


where 


N s = 


Wt 

W...C 


Equation (18) was derived from the simplified expression (eq. (9)) of the effectiveness of 
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a counterflow heat exchanger. The derivation of the lengths and terminal temperatures 
of the various regions is essentially completed. Before the pressure-drop function of the 
boiler can be expressed in closed form, however, the expression for the boiling -region 
pressure drop (eq. (7)) must be integrated. 


Integration of the Two-Phase Pressure-Drop Equation 

An expression for the local quality as a function of the length y is required to inte- 
grate equation (7). The amount of vapor generated for any length y, measured from the 
start of boiling, can be obtained by an energy balance on the heating fluid. The amount of 
vapor generated is proportional to the temperature drop in the heating fluid as given in 
equation (19): 


w v (y) 


w h°h 

H w 


< T h ' T hb) 


(19) 


The local temperature difference can be expressed in the form of the effectiveness 
relation of equation (12), and after appropriate substitution and simplification of terms 
the expression for w y (y) becomes 


w.c. 

w v (y) = JL5. (t. 


H. 


hs 


sat- 1 


exp (-N bLb j 


/N K y> 

I exp | -if- 1 - 11 


w 



( 20 ) 


or equivalently, 


x 



exp (-N bLb ) 



(21) 


The variable of integration in equation (21) is y; the length L b , as given by equation (14), 
is a constant and corresponds to the upper limit of the integral of equation (7). 

The integration of equation (7) can now be performed. The integral is divided into 
two parts corresponding to the plug and the swirl -wire or open -tube boiling regions as 
expressed in equation (22a): 
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(22a) 


f w 2 g rV"*P 
AP _ pl,b w w g p 

a b 

p w,v 

T 0 


/ 

Jo 


u 


x dy + 


Xk f‘b 

P ».v "V* p 


x dy 


Substituting the expression of equation (21) for x into this equation allows the closed- 
form integration of equation (22a): 


- V] 

b b \ *w, v 


N b< L pl - V - l } + H V* 


W, V 


exp [N b (L pl - L p) ] * fy L pl - L p )}j 


- Vb 


(22b) 


Rewriting the integrated expression in terms of dimensionless parameters results in 


Pw,v N b 


where 


(i • ‘) “(r^) * (i ‘ ‘) ■ ‘]j 


(23) 


r = y> bg p 

iw g 

b,^ & s 


V ■ N b< L P i - V 


and 
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<b° T IhS ' T Ttb =1 ' BiP K 1 ' 1 -) = 

^hs " 1 sat x ' 


w H X 
w w o 


w h c h< T hs - T sat> 


Equation (23) expresses the pressure drop of the boiling region as the frictional pressure- 
drop gradient for all gas flow in the swirl -wire or open-tube region, multiplied by 
X o yN b and the expression in braces. This latter expression can be considered as the 
dimensionless length weighting factor which incorporates the effects of the nonuniform 
heat-flux distribution. In addition, if r is 1 (i. e. , if the tube has no plug insert and 
a uniform internal geometry), the latter term in the braces is 0. Hence, this term is the 
increment of the weighting factor caused by the plug insert. 

This equation completes the derivation of the constituents of the boiler pressure 
drop. Next, these components are combined, and the boiler equations are summarized. 


Summary and Synthesis of Boiler Equations 

The total pressure drop of the boiler tube is the sum of equations (2), (3), and (23), 
as expressed in the following equation. 



Equation (24) can be expressed in dimensionless form by defining the pressure gradient 
in the superheat region at design flow times the total tube length as the normalizing fac- 

tor AP Vjd : 


AP. 


v,d 


_ g S f s w w,d*t 


v ,d 


(25) 


The resultant pressure-drop function becomes 
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( 26 ) 
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AP. 
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where the dimensionless parameters are 


W w= — 


w. 


w, d 
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p w, v 
p v, d 
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p,s 


^p g p p v, d 

f s g s p w,Z 


b,s 


_ f b, a 


The intermediate variables entering into the pressure-drop function can be expressed in 
the form of known or computable input variables as follows: 


L s = 1 - S - L p 


L s = 0 


V 1 ^ 1 

L P + ^ ^ 1 


(27) 


Lb= ^ln ( I -,,)- 1 


(28) 
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(29) 


e b = 


w w H w X o 


w h°h^ T hs " T sat ) 
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P 


d - V N P 


In 




(30) 


where 


and 


T - T 
_ sat wi 

P T t 
A hb ‘ wi 


a 


P 


w w c w, l 
w h®h 


N pl = N b( L pl " V 


L P1 >L P 


N pl = ° 


L, < L 
pi p 


(31) 


As mentioned previously, X Q in equation (29) is 1 except if either > 1 - L p or 

e b > 1. 0. In these cases X Q is adjusted so that the constraint = 1 - L p is satisfied. 

The temperature relations required for calculating the effectivenesses are expressed 
by equations (32) to (35): 


^hs ^hi “ a s^wo “ ^sat^ 


(32) 


T = T 
A hb A hs 


w H X 
w w o 

w h c h 


(33) 


T ho T hb " a p^ T sat " T wi^ 


(34) 
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T = T, . - (T. . - T 0 .\(l - ajexp -N (1 - ajL 
wo hi ' hi sat)' s' s' s' sj 


( 35 ) 


The saturation temperature is related directly to the saturation pressure of the boiler 
P ga f The saturation relation between T ga ^ and P ga ^ is expressed as 


^sat 


A 

B - U(P^ t ) 


(36) 


The parameters A and B were adjusted for each working fluid to fit its standard satur- 
ation curve. The saturation pressure, P gat was assumed to be that at the pinch point 
and is evaluated by equation (37): 

p sat = P o + ap b + AP s < 37 > 

The boiler exit pressure P Q depends on the exit boundary condition. For studying the 
potential for flow maldistribution among parallel tubes connected to a plenum, the most 
appropriate exit boundary condition is a constant exit-plenum pressure, which is indepen- 
dent of the flow variations in any one tube. On the other hand, for evaluating the pres- 
sure drop for uniform flow distribution in the tube bundle or for comparison with 
pressure-drop data from a boiler operating in a power loop with a turbine, the appropri- 
ate exit boundary condition is a nozzle-type restriction. A choked nozzle characteristic, 
as given by equation (38), was assumed for the latter situation: 

P = K, w X %[t (38) 

o tn w o f wo ' ’ 

For exit qualities less than 1, it was assumed that the nozzle pressure was determined 
by the gas flow rate w w X Q alone. 

The variation of the pressure level in the boiler may cause a significant variation in 
the vapor density. The variation of the dimensionless density factor R d was, for sim- 
plicity, related to the saturation pressure and temperature by the perfect gas law as 
expressed in 


R _ w, v 
K d ~ 


sat 


R_T 


p v,d "Vsat^r 


(39) 


where 
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p = Psat , d 

T sat, d^g 

was chosen so that is 1 at the reference design point. 

Conceptually, equations (24) to (39) constitute a complete mathematical model of the 
boiler operation. Because of the nonlinearity and complexity of the equations, a closed- 
form solution of the pressure drop in terms of the primary input variables T^, w , 

T hi> and cannot be obtained. Some insight into the interactions of the various equa- 
tions can be gained by representing them in block diagram form, as shown in figure 4. 

The boundary conditions are a specified working-fluid flow at inlet and a choked nozzle at 
exit. Superheat was assumed to be present at the exit. In the upper left of the figure, 
the pressure-drop equations are represented. The primary input variables to this section 
are the preheat length Lp and the effectiveness of the boiling region. These vari- 
ables are generated by the boiler heat-transfer equations as represented in the lower 
part of the figure. The major independent inputs to this section are the flow rates W w 
and and the inlet temperatures T w ^ and T^. The major dependent input to these 
equations is the saturation temperature T ga £ which is specified by the equations of the 
exit boundary condition and the pressure-drop relations as represented in the upper 
right of the figure. The major inputs to this latter section are the working -fluid flow, the 
pressure drop in the boiling and superheat regions, and the nozzle boundary condition. 

The interaction of the saturation pressure with the pressure drop, the staturation temper- 
ature, the density, and the outlet pressure is indicated by the various loops shown in 
figure 4. 

This discussion completes the derivation of the boiler equations and the defintion of 
their interaction. Before proceeding with the analysis of the theoretical pressure-drop 
function and its dependence on the various input variables and parameters, it is desirable 
to estimate the accuracy of the theory. This estimation was made by comparison with 
experimental data from a single-tube water boiler. 


COMPARISON OF THEORY WITH WATER-BOILER DATA 

The data presented herein were obtained under NASA contract NAS3-9420 from a 
counterflow boiler with water as both the heating and working fluids. This boiler is de- 
scribed in detail in appendix C. Tests were conducted around the design point (listed in 
table I) to establish the variation of boiler pressure, temperature, and pressure drop in 
response to variation in the working- and heating -fluid inlet temperatures and flows. The 
range of the variations was wide enough to determinethe dependence of the pressure drop 


16 



on the input variables and, hence, to permit a good comparison with theory. 

The theoretically predicted variations of the boiler pressure drop were calculated by 
using the equations derived in the previous section, with the parameters adjusted for the 
water boiler. The equations were programmed on a digital computer. A listing of the 
FORTRAN IV program and the calculation of the numerical values of the water -boiler 
parameters are given in appendix B. 

The comparison of the measured and predicted pressure drops as a function of the 
working -fluid flow is shown in figure 5(a). The pressure drops are shown in the dimen- 
sionless form normalized by the constant AP v ^ = 9.7 psia (67 000 N/m ). The flow 
rate is shown as a fraction of the design flow. Note the excellent agreement between the 
predicted and the measured data in both the level and variation of the pressure drop for 
flow rates below design. For flows higher than design, the overall variation of pressure 
drop was satisfactorily predicted. The largest error was in the transition from super- 
heated vapor to wet vapor, which occurs at the exit of the boiler as shown in figure 5. In 
the actual boiler, the transition is much less sharp than for the analytic model, because 
of the coexistence of droplets and superheated vapor at the exit. This nonequilibrium 
condition was encountered in the boiler testing as indicated by (1) high and low readings 
of thermocouples in the superheat region and (2) the decrease in superheat temperature 
that occurred in an adiabatic, coiled-tube section connected at the exit of the boiler. 

Under these conditions it could be expected that the predicted variation in exit superheat 
and quality would differ somewhat from that measured. However, the trend of the vari- 
ation is predicted reasonably well, as shown in figures 5(c) and (d). Shown in figure 5(b) 
is the expected and actual variation in the boiler -outlet pressure as a function of the flow 
with the assumption of a choked nozzle at the exit. In the transition region, the pressure 
against flow characteristic deviates from the linear relation expected for adiabatic gas 
flow. 

A further check on the validity of the pressure-drop analysis is the ability to predict 
the variation of pressure drop caused by changes in the heating -fluid flow and inlet tem- 
perature, as well as the subcooled inlet temperature of the working fluid. The agreement 
of the pressure drops both in level and variation is excellent, as shown in figure 6. The 
predicted values generally fall about 10 percent below the measured values. This dif- 
ference was expected because of the neglected effects in the simplified homogeneous two- 
phase flow model used in the pressure-drop prediction. 

This comparison of theory and data shows that the analytic model accurately predicts 
both the pressure-drop level and its variations in the single-tube boiler. The analytic 
model should therefore be useful in the analysis of the pressure -drop characteristics that 
cause the maldistribution in the multitube configuration. Although the theory is shown to 
predict boiler performance for only the particular data available, the analysis might be 
valid over a much broader range. Applicability to other ranges of boiler operation and 
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design points (such as boilers designed for low pressure drop or high heating -fluid flow 
rates) has not been evaluated. There appears to be a paucity of data in the literature for 
total performance. 

EXAMINATION OF BOILER PRESSURE-DROP CHARACTERISTICS 

The analysis was performed to establish the general characteristics of the 
counterflow -boiler pressure-drop function and how the various parameters and input vari- 
ables influence this characteristic. In view of the complexity of the function and the 20 
or more variables involved, the analysis was limited to the effects of the parameters 
r s,J> V’ and U b and the input variables T wi , w w , w fa , and T^. 

To make the analysis more applicable to present-day space power systems, the ref- 
erence values for the boiler parameters were chosen to match those of a seven-tube mer- 
cury boiler presently being constructed for the SNAP-8 system. This boiler is predicted 
to have a multivalued pressure-drop function which is to be made single valued by inser- 
tion of a nozzle at the entrance of each tube. The influence of the nozzle on the pressure- 
drop characteristic is treated in the analysis. A brief description of the mercury boiler 
and how it differs from the water boiler is given in appendix D. In this analysis both 
choked-nozzle and constant-pressure boiler -tube discharge conditions are considered. 

The pressure-drop function for a choked-nozzle discharge appropriate when considering 
the boiler as a whole is computed by using the FORTRAN IV program of appendix B with 
constants adjusted for the mercury boiler. The pressure-drop function for constant exit 
pressure appropriate when considering one tube acting in a multitube configuration is 
computed by using the FORTRAN IV program given in appendix E. 


Components of Pressure-Drop Function 

Figure 7(a) illustrates the typical shape and composition of the multitube boiler 
pressure-drop-against-flow function. The boundary condition at the exit is a choked 
nozzle. The total dimensionless pressure drop <p^ (solid curve) is the sum of the com- 
ponent pressure drops (dashed curves) of the preheat, boiling, and superheat regions. 
The pressure drop due to an inlet orifice is not included. For all but the highest flows, 
the preheat pressure drop is negligibly small so that only the boiling -region and 
superheat-region pressure-drop components must be considered. The latter exhibit a 
complex variation with flow, which can be explained by considering the variation of the 
lengths of the regions in relation to the plug and tube lengths. These lengths are shown 
in figure 7(b). At very low flow rates, the preheat and boiling regions are so short that 
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substantial superheating occurs in the plug with the attendant high velocities and a gas 
pressure drop in the plug. This pressure drop at first increases with flow. However, at 
the same time, the preheat and boiling lengths increase to gradually supplant the super- 
heat region in the plug and, thus, produce a net decrease in the superheat pressure drop. 
The boiling pressure drop, in the meantime, increases rapidly so that the total pressure 
drop continues to rise until a flow of 0. 6 is reached. At this point the superheat region 
has been moved completely out of the plug, and the boiling pressure drop has reached a 
maximum. As the flow rate is increased further, the high-quality end of the boiling 
region is pushed out of the plug region, which causes a definite decrease in both the boil- 
ing and total pressure drops. This negative 3 <p^/ 3 W w persists until approximately de- 
sign flow for the case shown. The superheat pressure drop in the meantime is increasing 
gradually as a result of the increasing flow but is overcome by a decreasing superheat 
length as the boiling and preheat regions continue to expand. (In other words, above 
design flow the superheat length begins to shrink so rapidly that a decrease in superheat 
pressure drop occurs despite the increase in flow. ) The total pressure-drop function 
shows a corresponding increase in the negative slope at these flow rates. As the boiler 
runs out of superheat length, the presssure drop reaches a minimum. Any further in- 
crease in flow causes a gradual increase in the boiler pressure drop. This increase is, 
however, not very rapid because of the decreasing exit quality and length of the boiling 
region. If the preheat region expands rapidly enough, a second region of negative slope 
could result at this point. This region would correspond to the situation widely treated 
in the literature (refs. 2 to 4 and 9). This second region of negative slope is, however, 
not likely to occur in space power boilers because of the greater liquid-pressure-drop 
gradient in the plug region. Furthermore, in these boilers, inlet orifices are usually 
employed to add an overriding positive slope above design flow. 

As shown in the preceding analysis, the plug insert is primarily responsible for the 
negative 3 <p^/ 3W w characteristic. Removal of the plug does not necessarily solve the 
problem, however, as can be seen in figure 7(c). Here the r parameter has been 

Sy 

set equal to 1, which effectively eliminates the plug insert (see eq. (19)). Yet the 
pressure-drop function exhibits a negative -slope region, even though this slope as well as 
the pressure-drop level is much reduced. Thus, two distinct mechanisms are at work. 
One is the result of the boiling and superheat regions moving out of the plug insert. The 
other occurs when the superheat length is rapidly eliminated from the tube as the boiling 
region expands with flow. 


Effect of Boiling Pressure Level and Its Dependence on the Exit Boundary Condition 

The simple explanation given in the preceding section is not completely adequate in 
describing the magnitude of the negative dcp/dW^. This slope is substantially influenced 
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by the variation of the vapor density as the flow rate is changed. This variation is caused 
by the choked nozzle assumed at the exit of the boiler tube bundle which results in a linear 
increase in exit pressure with flow as long as superheat remains at the exit. The higher 
pressure causes a higher vapor density, which acts to reduce the pressure drop (see den- 
sity parameter in eq. (26)). In contrast, if the pressure in the boiler were fixed, 
the d<p/d W w would be less negative. Such a situation is actually of primary interest in 
evaluating the possibility of a single tube operating a different flow but at the same pres- 
sure level and pressure drop as other tubes in a bundle. Such a tube is best described by 
a constant exit pressure corresponding to the exit-plunum pressure as determined by the 
total flow through the boiler. Such a boundary condition describes only the simplest pos- 
sible maldistribution in a tube bundle, which in general can take on a multiplicity of flow- 
distribution modes. Both the simplest mode and these more complex modes can exist 
if the single-tube pressure-drop-against-flow function is multivalued. Hence, attention 
was focused on the pressure-drop characteristic of a simple single tube with constant 
exit pressure to evaluate the actual potential for maldistribution. The characteristic is 
evaluated with the assumption that the heating fluid in the tube bundle is unmixed (i. e. , 
each tube has its individual supply of heating fluid). This assumption represents the 
worst case, in that the slope dtp/ dW is more negative for the case of unmixed flow than 
for mixed flow, where the variations of heating -fluid temperature in response to working- 
fluid flow changes are minimized. 

The pressure-drop function <p^ as a function of W w for various pressure levels is 
shown in figure 8 (solid curves). The tube bundle with a choked-nozzle characteristic 
(dashed curve) was added for contrast. The difference between the dashed and the solid 
curves, as well as the spacing between the solid curves, indicates the substantial effect 
of the pressure level on density. The general multivalued characteristic is evident at all 
pressure levels, but as expected, the slopes of the curves are less than the choked- 
nozzle characteristic shown. Of particular interest is the single-tube pressure-drop as 
a function of flow around the design condition (P Q /P o ^ = 1. 0). By coincidence, the 
pressure-drop function for this pressure level is at a maximum at design flow (W = 1. 0). 
This operating point is unstable for the single tube operating in a tube bundle (i. e. , con- 
stant pressure at exit) and would probably lead to wide maldistribution of flow. Under 
such conditions, a mixture of superheated vapor and liquid would emerge from the boiler 
exit plenum. Similar results would be obtained if operation at other tube-bundle pressure- 
drop levels were to be attempted, as evidenced by the multivalued nature of the single - 
tube functions at all exit pressure levels. 
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Eliminating Negative-Slope Region in Pressure-Drop Function 


The most direct way of eliminating the negative -slope region in the pressure-drop- 
against-flow function is to increase the liquid flow resistance. This resistance can be 
increased by narrowing the flow passages or by increasing the spiral pitch in the plug 
(i.e., increasing the r g ^ parameter). This procedure, however, also increases the 
resistance to flow in the plug boiling region (increases the r parameter) which negates 
any improvement resulting from the increased liquid flow resistance. Designing a tapered 
flow passage or using a plug with a tight pitch in the preheat region and loose pitch in the 
boiling section has been attempted, but these more complicated cases will not be treated 
herein. A much simpler approach, which is to be used in the latest design of the SNAP- 8 
boiler, is to insert an orifice at the entrance of ea,ch tube. The only limitations to this 
method are the allowable pumping loss and the possibility of plugging if the orifice is 
made too small. 

The effect of the inlet orifice on the pressure-drop characteristic is illustrated in 
figure 9. The dimensionless parameter r^ is the ratio of the orifice pressure-drop 
parameter AP y The solid curves represent the characteristics of a tube with con- 
stant exit pressure (i. e. , a single tube operating in a large bundle); the dashed curves are 
for the boiler with a choked nozzle at the exit. The solid curves show that an orifice of 
r^ = 0. 5 is required to prevent a multivalued pressure drop for the single tube (constant 
exit pressure). This value is therefore the minimum required to eliminate the potential 
for maldistribution in a tube bundle. An orifice of greater resistance (r^ > 1. 0) is 
required to eliminate the negative dcp^/ 9W w region for the multitube boiler with a choked 
nozzle (dashed curves). In general, the amount of orif icing is dependent not only on the 
boundary condition, but also on the plug geometry and heat flux as well. 


Effect of Plug Geometry 

The plug length L j and its flow resistance per unit length, as defined by the r g A 
parameter, have a substantial effect on the pressure-drop function. This effect is illus- 
trated in figure 10. The pressure drops shown include the component due to an orifice of 

r. = 0. 5. As the r parameter increases (fig. 10(a)), the pressure drop and the peak 
i ** 

in the pressure-drop curve increase. Thus, although an inlet orifice of r^ = 0. 5 was 
sufficient to eliminate the maldistribution for the boiler with r_ = 36. 0, a boiler with 
r c = 100. 0 would require an inlet orifice of r. = 0. 8. The dashed curve in figure 10 
represents the part of the pressure drop caused by an orifice with r^ = 0. 5; higher 
values of r^ produce proportionately higher pressure drops and steeper positive slopes. 
The effect of the plug length on the shape of the pressure-drop function is illustrated in 


21 


figure 10(b). The increase in pressure drop is substantial (about 40 percent at design 
flow) when the plug length is increased from 15 to 25 percent of the tube length. As with 
the r parameter, an increase in the plug length requires an increase in the size of 
the liquid-pressure-drop orifice to eliminate the negative -slope region. Conversely, 
shortening the plug tends to reduce the amount of required orificing. It is not desirable, 
however, to remove the plug completely because of heat-transfer considerations. 


Effect of Heat-Flux Distribution on Pressure-Drop Function 

It has been demonstrated that the density and plug-geometry factors affect the shape 
and magnitude of the pressure-drop function. A third category of parameters, associated 
with the heat-flux distribution and rate, must yet be considered. It has already been 
demonstrated in the section comparison of theory with water-boiler data, that the pres- 
sure drop varies with the input variables that influence only the heat transfer (i. e. , the 
heater-loop flow and inlet temperature, see figures 6(a) and (b). In addition to changing 
the pressure-drop level, these variables also change the shape of the <p t (W w ) function, 
as shown in figures 11(a) and (b). 

Figure 11(a) shows that as the heating -fluid flow increases the pressure drop in 

the boiler increases at all power loop flow rates. This increase results from the 
increase rate of vaporization and shorter boiling and preheat lengths associated with the 
larger average temperature difference. The variation of length is indicated indirectly by 
the shift in the minimum point of the pressure-drop curve as increases. This min- 
imum point indicates the working -fluid flow for which the superheat length becomes zero. 
As expected, the working -fluid flow at which the boiler runs out of superheat length 
increases with increasing heating -fluid flow. This increase has a desirable effect on the 
shape of the pressure-drop function in that it extends the region of positive slope. At a 
heating flow of twice design (W^ = 2), the negative slope emerges only at flow rates 
above 1. 6 of design working -fluid flow. 

A much smaller liquid-flow resistance (larger orifice) would be needed to eliminate 
this negative slope than is required at design heating -fluid flow (W^ = 1. 0). Thus, an 
increase in heating -fluid flow, if the additional pumping power can be justified, will 
generally produce an improvement in the pressure-drop characteristic. In the limit as 
W h becomes large, the heating -fluid temperature drop becomes negligibly small; hence, 
the curve = 100. 0 represents a boiler with practically uniform heating-fluid temper- 
ature. For this case, the region of steep negative slope in the superheat component of 
the pressure drop (see fig. 7(b)) represents a high working -fluid flow (relative to design). 
This causes the liquid-pressure-drop resistance to become large enough, even without 
an orifice, to prevent the total pressure-drop function from having a region of negative 
9<?/3W w . 
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The effect of the heating -fluid inlet temperature, as shown in figure 11(b), is basi- 
cally the same as that of the working-fluid flow. The temperature variation is expressed 
in the form of a dimensionless parameter (T^ - T ga ^. ^)/ (T^ ^ - T gat ^). As T^ 
increases, the average temperature difference between the heating fluid and the saturation 
temperature T ga j. increases and produces a greater pressure drop and a more positive 
9^/3 W w at W w - 1.0. 

The third factor which can change the heat flux and hence the pressure-drop char- 
acteristic is the variation of heat-transfer coefficients. Figure 11(c) shows the effect of 
variation in the boiling coefficient. Note the curious effect that better heat-transfer 
produces a less desirable characteristic as far as potential maldistribution is concerned. 
In contrast, for a U b of one-fifth of design (U b = 100), the rate of vaporization is so 
slow that no peak in the pressure-drop function occurs. Such a boiler, however, would 
generate vapor without superheat at a quality of 60 percent at design liquid flow rate. 

Thus, decreasing the heat-transfer coefficient does not appear to be a good method of 
eliminating the undesirable pressure -drop characteristic. 

A question arises at this point as to why increasing the heat flux by increasing 
heating -fluid flow and inlet temperature helps to reduce the negative -slope region in the 
pressure-drop curve, while the increase in heat flux caused by an increased boiling coef- 
ficient produces the opposite effect. The basic difference is apparent in the equation of 
the boiling length (eq. (14)) as expressed in terms of the primary variables: 


1 

, w w H w X o 

_ "h^hs - T sat) _ 

Increasing shortens the boiling length but does not affect the effectiveness term inside 
the parentheses. In contrast, increasing heating -fluid flow or temperature decreases 
the effectiveness. (Note that T h s ~ T hi' ) 1° the other direction, as the effectiveness 
approaches 1, the boiling length becomes large and very sensitive to working -fluid flow. 
This situation leads to a very rapid expansion of the boiling length and a corresponding 
decrease of the superheat length as the working -fluid flow increases. It is this effect 
that is responsible for the negative slope of the pressure-drop function of a tube without 
the plug insert. 

The three variables W h , T hi , and U b are, of course, not the only ones that influ- 
ence the pressure-drop function through changes in heat flux. They are, however, the 
major ones. The subcooled- working -fluid inlet temperature T w ^ has little effect on the 
pressure-drop function, both for the water boiler (see fig. 6(c)) and the mercury boiler 
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(not shown here but similar). This condition is primarily caused by the high values 
(see table n) which make the preheat length insensitive to changes in the inlet temperature 
(see eq. (10)). For smaller Np values (e.g., for low values of Up), the pressure drop 
would become more sensitive to the amount of subcooling. This situation is not typical, 
however, in liquid-metal space power boilers, particularly those that have plug inserts to 
improve the heat transfer. 


CONCLUDING REMARKS 

A theoretical boiler model was developed and verified by comparison with experi- 
mental data. An analysis based on this model has led to a number of conclusions regard- 
ing (1) the factors involved in producing the multivalued function of pressure drop against 
working -fluid flow for a boiler tube and (2) the possible approaches to boiler design that 
would make the multivalued function single valued and thus eliminate this potential for 
flow maldistribution in a tube bundle. These design approaches are directed toward 
eliminating the negative -slope region of the function (i. e. , the region where pressure 
drop decreases with increasing flow). 

The analysis has indicated that this negative -slope region is the result of a combina- 
tion of performance factors. In a tube with a plug insert, the primary factor is the rapid 
decrease in the boiling pressure drop in the plug region as the working -fluid flow in- 
creases. Another factor, which is present in a boiling tube even without a plug insert, is 
the decline of the superheat length and pressure drop as the boiling region expands with 
increasing flow rate. Still another factor is the variation of vapor density with changes in 
boiler pressure level, where the pressure level is determined by the boiler -exit boundary 
conditions. The boundary condition for a single tube operating in a large tube bundle is 
that exit pressure be nearly independent of the flow in the individual tube (i. e. , constant 
pressure). For this boundary condition, the pressure drop in the single tube decreases 
less with increased flow than for a boundary condition where exit pressure increases 
directly with flow (i. e. , choked nozzle). 

The analysis investigated the effect on the negative -slope region of the variation in 
some of the geometric and thermodynamic parameters of the boiler. These parameters 
show varying degrees of effectiveness in eliminating the negative-slope region. Mini- 
mizing the plug length or its resistance reduces the peak in the pressure -drop function 
as well as the steepness of the negative-slope region. Such a reduction of the plug is 
limited, however, by heat-transfer considerations. Increasing the heating-fluid flow or 
inlet temperature, and thereby the average temperature difference, also helps to make 
the slope of pressure drop against working-fluid flow less negative. In contrast, in- 
creasing the heat-transfer coefficient in the boiling region increases the maximum in the 
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pressure-drop function and also increases the negative slope. Adding an orifice at the 
entrance of each tube is a simple method of assuring that pressure drop always increases 
with increasing flow. The amount of orificing required increases with the steepness of 
the negative slope. 

All the proposed methods of improving the pressure-drop characteristic, except the 
plug modification, require increased pumping power or higher heating temperatures and 
thus result in a system performance penalty. The reduction of the plug length or plug- 
flow resistance to the minimum consistent with good heat-transfer performance appears 
to be the best method of improvement. Nevertheless, even with the plug completely 
removed, a negative slope in the pressure-drop function may still exist and hence require 
an orifice to eliminate the potential for flow maldistribution. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 7, 1967, 

701-04-00-02-22. 



APPENDIX A 


SYMBOLS 


Mathematical 

FORTRAN 

Description 

symbol 

variable 


A 

AGAS 

constant in eq. (36) adjusted to fit saturation- 
temperature-against -pressure curve for working- 
fluid cross-sectional area of boiler tube 

A 

c 


cross-sectional area of boiler tube 

A pl 


cross-sectional area of plug flow passage (perpen- 
dicular to flow), ft 2 ; m 2 

A s 


cross-sectional area of tube in superheat region, 
ft 2 ; m 2 

a p 

APH 

ratio of working -fluid to heating -fluid heat capacity 
in preheat region, w w c w , 

a s 

ASH 

ratio of working -fluid to heating -fluid heat capacity 
in superheated region, w w c w v / w h c h 

B 

BGAS 

constant in eq. (36) adjusted to fit saturation- 
temperature-against-pressure curve for working 
fluid 

c h 

CH 

specific heat of heating fluid, Btu/ (lb)(°F); J/ (kg)(°K) 

c w l 

CPP 

specific heat of working -fluid liquid, Btu/ (lb)(°F); 



j/(kg)(°K) 

V 

CPV 

specific heat of working -fluid vapor, Btu/ (lb)(°F); 

W, V 


J/(kg)(°K) 

D 


hydraulic diameter, ft; m 

D. 

1 


inside diameter of boiling tube, ft; m 

D o 

DBL, DPH 

outside diameter of boiling tube used for heat- 
transfer calculations, ft; m 

f 


friction factor 

f pl,b 


friction factors in plug boiling region 

f b,d 

26 


friction factors in swirl -wire region 
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Mathematical 

FORTRAN 

Description 

symbol 

variable 


f p’ f b’ f s 


friction factors in preheat, boiling, and superheat 
regions 

f pl,P 


friction factors in plug preheat region 

f pl,v 


friction factors in plug for all vapor flow 

G p’ G b’ G t 

RPHDP, RBLDP, 

dimensionless equivalent length factors for calcu- 


and RTDP 

lating the preheat, boiling and, total pressure 
drops 

g 


acceleration of gravity, ft/sec^; m/sec^ 

gp> g s 


geometry factors in pressure-drop equations of pre- 
heat, boiling, and superheat regions, respectively 

H w 

HV 

heat of vaporization of working fluid, Btu/lb; J/kg 

*Sv, s 


heat-transfer coefficient of working fluid in super- 
heat region 

K tn 

TNK 

choked turbine nozzle constant at boiler exit, 
^*od^ w w, d^o V T wo, d 

l 


generalized length along tube, ft; m 

*b’ L b 

EBOIL 

boiling length, ft; m (dimensionless, 

Z ,L 
P’ P 

ELPH 

preheat length, ft; m (dimensionless, 

W 

EPLUG 

plug length, ft; m (dimensionless, l 1 ^.) 

^s ,L s 

ELSH 

superheat length, ft; m (dimensionless, 1-/1*.) 

*t> L t 

ALT 

total boiler tube length, ft; m 

N 


number of 90° turns for spiral plug 

NTU 


number of transfer units 

N b 

BNTU 

NTU of boiler based on U. 

b 

N P 

PHNTU 

NTU of boiler based on U 

p 

V 

PLNTU 

NTU based on boiling region in plug, N.^(L j - L^) 

N s 


NTU of boiler based on U 

s 

n 


number of turns per inch 

P 


pressure, psia; N/m abs 
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Mathematical 

symbol 


FORTRAN 

variable 


Description 


p i 

PIN 

P o 

POUT 

P o,d 


P sat 

PSAT 

P sat, d 



PSA TO 


AP p> AP b> AP s 


A p t 

DP 

AP v,d 

AVPRS 

AP. 

1 


P 


P P1 


p s 


R d 

ROW 

Re s 


Re pl 


R g 

RGAS 

r b,s 

RBSP 

r i 

PDI 

r P,s 

RPH 


n 

working-fluid inlet pressure, psia; N/m abs 

working -fluid outlet pressure, psia; N/m abs 

working -fluid outlet pressure (design value), psia; 
N/m ^ abs 

working-fluid saturation pressure, psia; N/m abs 
working -fluid saturation pressure (design value), 

o 

psia; N/m abs 

working -fluid saturation pressure used for initial 

p 

estimate in iteration, psia; N/m abs 
pressure drops in preheat, boiling, and superheat 

o 

regions, psi; N/m 

2 

total boiler-pressure drop, psi; N/m 

normalizing pressure-drop factor based on pressure 
gradient in superheat region at design flow, psi; 
N/m 2 

orifice pressure drop at tube inlet at design flow, 
psi; N/m 2 

wetted perimeter, ft; m 

wetted perimeter in plug flow passage, ft; m 

wetted perimeter in superheat open-tube region, 
ft; m 

dimensionless vapor density, p w y / P r 
Reynolds number in superheat region 
Reynolds number in plug region (liquid) 
working -fluid gas constant 
ratio of boiling to superheat friction factors 
ratio of orifice pressure drop to AP , 

V , u. 

ratio of preheat-plug-region pressure gradient to 
superheat pressure gradient 
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Mathematical 

symbol 

FORTRAN 

variable 

Description 

r s,u 

RPS 

ratio of pressure gradients for vapor at rated flow in 
plug region to those of vapor at rated flow in swirl - 
wire region 

T hb 


heating -fluid temperature at point where boiling 
starts, °F; °K 

T hi 

THL 

heating -fluid inlet temperature, °F; °K 

T hi, d 


heating -fluid inlet temperature (design value), °F; 
°K 

T ho 


heating -fluid outlet temperature, °F; °K 

T hs 

THLSH 

heating -fluid temperature at point where superheat- 
ing starts, °F ; °K 

T sat 

TSAT 

working-fluid saturation temperature 

T sat, d 


working-fluid saturation temperature (design value) 

T . 

W1 

TPO 

working-fluid inlet temperature, °F; °K 

^wi, d 


working -fluid inlet temperature (design value), °F- 
°K 

T 

wo 

TPL 

working -fluid outlet temperature, °F; °K 

T ws 


working -fluid saturation temperature at start of 
superheating, °F; °K 

u b 

UBL 

overall average heat-transfer coefficient in boiling 
region, Btu/(ft 2 )(hr)(°F); j/(m 2 )(hr)(°K) 

u h,d 


overall average heat-transfer coefficient in boiling 
region (design value), Btu/(ft 2 )(hr)(°F); 
J/(m 2 )(hr)(°K) 

U P 

UPH 

overall average heat-transfer coefficient in preheat 
region, Btu/(ft 2 )(hr)(°F); J/(m 2 )(hr)(°K) 

U s 

USH 

overall average heat-transfer coefficient in super- 
heat region 

V 


velocity of working -fluid flow, ft/hr; m/hr 

w h’ W h 

WH 

heating -fluid flow (dimensional, lb/hr (kg/hr); 


dimensionless, = w j 1 / w } 1 ^ 
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Mathematical 

symbol 

FORTRAN 

variable 

Description 

w h,d’ W h,d 

WHDG 

heating -fluid design flow (dimensional, lb/hr 
(kg/hr); dimensionless) 

w v (y) 


local vapor flow along boiling length y, lb/hr; kg/ hr 

w w’ W w 

WP 

working -fluid flow (dimensional, lb/hr (kg/hr); 
dimensionless), W w = w^w w d 

w w,d’ W w,d 

WPDG 

working -fluid design flow (dimensional, lb/hr 
(kg/hr); dimensionless) 

X o 

QUAL 

outlet quality at end of boiling region 

X 


local quality along boiling length y, lb/hr; kg/hr 

y 


length along boiling region, ft; m 

a(AP)/aw w 


slope of boiler-pressure-drop-against-working -fluid- 
flow function 

V 3W w 


slope of dimensionless pressure-drop function 

e b 

BBOIL 

effectiveness of boiling region 

€ P 

BPHT 

effectiveness of preheat region 

e s 


effectiveness of superheat region 

M 


viscosity 

P 


density, lb/ft 3 ; kg/m^ 

Pr 

ROWRE 

working -fluid vapor density at design saturation 

3 3 

pressure used to normalize eq. (39), lb/ft ; kg/m 

p v,d 


working -fluid vapor density at design flow in super- 
heat region, lb/ft kg/m 3 

p w, l 

ROWL 

3 3 

working -fluid liquid density, lb/ft ; kg/m 

P W, V 

n 

RD 

O 

working-fluid vapor density, lb/ft ; kg/m 
boiler dimensionless pressure-drop function, 

a V iP v,d . 
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APPENDIX B 


DIGITAL PROGRAM FOR BOILER MODEL WITH CHOKED-NOZZLE DISCHARGE 
AND CALCULATION OF WATER-BOILER CONSTANTS 

Computer Program 


The equations of the analytical model were programmed on a digital computer and 
solved simultaneously by iteration. (See program listing at the end of this appendix. ) 

The iteration consists simply of adjusting the saturation pressure until the boiler-outlet 
pressure corresponds to the working -fluid flow under the conditions specified by the 
choked-nozzle boundary condition. A secondary iteration loop is provided to converge on 
outlet quality for the cases where the quality at the boiler exit is less than 1. The main 
body of the program is applicable to any counterflow boiler. However, numerical values 
and data statements in the listing apply specifically to the water boiler. The dimension- 
less parameters and the flows, pressures, and temperatures at design conditions re- 
quired for the evaluation of the boiler equations are listed in table I. The method of cal- 
culating the pressure-drop parameters (P 17 A , r r and r, ) is shown in a later 

Vj Cl Pj S Sj (* Dj S 

part of this appendix. The other parameters in the table can be calculated directly from 
the fluid properties and the specified design values of the input variables. The heat- 
transfer coefficients for evaluating N N, , and N were calculated by using the Dittus- 
Bolter equation for forced convection. The convective coefficient on the working -fluid 
side was based on the hydraulic diameter. The coefficient in the plug preheat region was 
multiplied by the factor (f pl p / f smooth tube^° 28 = 1 ‘ 19 > as developed from reference 12, 
to account for the improvement resulting from spiral flow. The convective coefficient in 


the superheat region was multiplied by a similar factor, (h w s /h smoot j 1 t^e) — 2.0, as 
obtained from figure 10 of reference 12 for a wire -coil -pitch to coil -diameter ratio of 5. 


The average boiling coefficient, based on previous experience with a water boiler, was 
known to be very high in comparison with the heating -fluid coefficient and the wall con- 


ductivity. For purposes of calculation, this boiling coefficient was assumed to be 6000 
Btu per hour per square foot per °F (34 000 J/(sec)(m^)(°K)). 


Calculation of Pressure-Drop Parameters for Single-Tube Water Boiler 

The fluid properties used for the calculations were taken at the design pressures and 
temperatures and were assumed to be constant. 
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Calculation of AP V d . - The normalizing factor for the pressure-drop function is 
defined as 


AP 


v,d 



which represents the pressure gradient in the superheat region at design flow multiplied 
by the total tube length. The friction factor f is based on the flow geometry of the 
superheat region, which includes not only the spiral wire attached to the inside tube wall, 
but also, in the center of the tube, a 0. 25-inch- (0. 625-cm-) diameter rod used to support 
the internal instrumentation. 

The flow cross-sectional area is 


A 


s 



£ 

4 


(0. 0462) 2 - (0. 208) 2 


0. 1944 in. 2 (1.25 cm 2 ) 


The wetted perimeter is 

p s = n(D i + D rod ) = ir(0. 0462 + 0. 0208) = 0. 21 ft (0. 064 m) 

The Reynolds number for the annular flow passage at design vapor flow is 

4w 

Re = d = 23 600 

s p u 0.210x0.0364 
s 


The smooth-tube friction factor at Re = 23 600 is approximately 0. 025. The friction- 
factor multiplier to account for the swirl wire is based on the data of reference 12. 

Figure 15 of this reference gives the friction-factor multiplier as a function of wire-coil- 
pitch to wire-diameter ratio s/D for various values of wire diameter. These data were 
extrapolated to the value of D w t re / D tube = 12 at S / D = 5 of the water-boiler tube, 

which gave a value of f g / f smoo tti tube ~ 17 * Now AP V d can be calculated. 


AP 


v , d 



<d' 


v, d 


t 
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AP 


v,d 


17. 0(0. 025) 


144 


lb/ft d _ 
psi 


0.21 


[8(32. 2)(0. 00135)'; 


(1J56x10 -4 )10. 0 


0. 155 


= 9. 7 psi (66 930N/m 2 ) 


Calculation of r - The definition of r is 
P»_i p, s 


r 


P, s 




The values of 
The factor d£ 
of 1. 


fp, A pl and p^ were computed on the basis of flow along the spiral path, 
pj/dl modifies the pressure gradient to correspond to the axial tube length 



dt cos 6 cos 75. 80 


4. 09 


The calculations are based on the geometry shown in the following sketch: 

Flow 



The wetted perimeter of the area normal to flow is 

Ppl = 2 [(sin e ) 0. 030 + 0. 00417] = 0. 0664 ft (0. 02024 m) 

The cross-sectional area normal to flow is 

A j = sin 75. 8X0. 00417x0. 030 = 0. 000121 ft 2 (10.89xl0 -6 m 2 ) 

The Reynolds number for design water flow is 
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4w. 


Re 


£,pl 


w, d 




4. 45 
0.0664(0.635) 


= 4270 


The smooth -tube friction factor at an Re of 4270, with the assumption of turbulent flow, 
was approximately 0. 038. The friction -factor multiplier which accounts for the bends in 
the spiral flow passage was calculated by use of a simplified form of the relation for 
equivalent length-diameter ratio for bends greater than 90° (ref. 13). (Symbols are the 
same as used in ref. 13. ) 




where D is the hydraulic diameter, R^ is the resistance due to the length of one 90° 
bend, R^ is the bend resistance due to one 90° bend, the total resistance R^ = R^ + R^, 
and N is the total number of 90° bends in a coil. The resistance due to the bending of 
flow can be extracted from the preceding expression. 


isa = (n + 1) 

D 2 


The number of 90° turns N for the spiral plug is 


N = 4nZpi = 258 

where N is the number of turns per inch (per cm). When this number is large enough so 
that (N + 1) = N, the following expression for equivalent length in terms of plug length 
l pi due only to the bending of flow can be written: 




4,1 w 


The total equivalent length of the plug can now be expressed as 


4nZ ,R,D 

l, = 4. 09 Z , + P 1 D 

pt pi 2 

/ 4nR, D 

= 4.09* j 1 + °- 

pl \ 4.09x2, 
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Given the hydraulic diameter D = 4A C /P S = 0. 0073 feet (0. 00226 m) and = 8 for 
r/D = 2. 88, the equivalent length becomes 


v =4 - o9 v r + 


120. 8x0. 0073) 
8. 18 I 


= (4.09 X 1.86)Z 


Pi 


All the numbers required for calculating r are now known: 

P> b 


r = 71 P L !p Py . d 
P,S M ‘s <• l p s 


M 

\v/ 


= 4 09 0. 038X1.86 0. 155 0.0664/0. 00135 \ _ Q ?8 
0.025x17 59.4 0.21 lo. 000121/ 


Calculation of r g a parameter. - The definition of r g 4 parameter is 

r - f Pl.b P p/M 

S ’' f b, S di p s\V 

The friction factor fp^ ^ for boiling in the plug is assumed to correspond to a single- 
phase gas Reynolds number of approximately 100 000 at design flow. The friction -factor 
multiplier is the same used for f p in the previous section (i. e. , 1. 86). The friction 
factor in the swirl-wire region was assumed to be equal to f_. All the other parameters 
are already known; hence, the value of r g a can be calculated: 


0.019x1.86 4 Q9 0.0664 /0.00135 


s ’ 4 0.025x17 


0.21 \0. 000121, 


= 150. 0 


The r fe g parameter defined as f^ Ji Q is 1 since f fe 4 was assumed to be equal to 


f s* 
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PROGRAM LISTING 


SIBFTC BLCN 


LIST 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


30 

31 

32 

33 

34 

35 

36 
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THIS IS A FORTRAN IV PROGRAM FOR COMPUTING THE PRESSURE DROP 
OF A SINGLE-PASS COUNTERFLOW BOILER. 

THE EXIT BOUNDARY CONDITION IS A CHOKED NOZZLE 

THE PROGRAM IS GENERAL EXCEPT FOR THE INITIALIZATION SECTION AND 
THE DATA CARDS. THE NUMERICAL CONSTANTS ARE SPECIFIED IN ENGLISH 
UNITS , ALTHOUGH ANY CONSISTENT UNITS CAN BE EMPLOYED. 

(ENGLISH UNITS ARE IN FT HR LB BTU PSI DEG-F ) . 

THE NUMERICAL VALUES ARE FOR THE SINGLE TUBE WATER BOILER 
THE DATA - AWP - SPECIFIES VALUES OF WORKING FLUID FLOW (IN UNITS 
OF FRACTION OF DESIGN FLOW ) FOR WHICH THE PRESSURE DROP FUNCTION 
IS TO BE EVALUATED. 

DIMENSION AWP (25) 

DATA AWP/ 1.0,0.02,. 06 *0.1, 0.2*0. 3*0.4, 0.5, 0.6, 0.7, 0.8,0. 9,1.1, 
11.2,1.3,1.4,1.6,1.8,2.0,0.95,1.05,1.15,1.25,1.35,1.5 / 

FORMAT (2I12?1F12.7,I12 ) 


WPO 


FORMAT ( 6F 1 2 • 7 ) 

FORMAT ( 49H 

125H THL PO 

FORMAT (20H1 INPUT PARAMETERS 
1 6H DK= , E 1 2 . 5 , 

FORMAT ( 1H ,9 ( E 1 2 . 5 ) ) 

FORMAT ( 50H1 WH/FNBL WP/FNPH 
1 60H UPH/TNK USH/AVPRS UBL/A1 
FORMAT ( 1H , 7 ( El 2 • 5 ) ) 

FORMAT ( 50H APH 

1 60H SHNTU PLNTU 

106 FORMAT (50H WP/PIN 
1 60H DP/ELPH TPL/ELSH 

107 FORMAT (1 H , 9IE12.5)) 

108 FORMAT (1H , 9(E12.5)) 


WH 


TP- 


PSATO 

/ 6F12.7 ) 

,6HICASE=, I5/7H KCMAX=,I5, 


ASH 

ROW 

QUAL /WH 


6H KPR = ,15 

RPH/FNSH 

ALT/A2 

PHNTU 


) 

RPS/RBSP 

EPLUG/A3 

BNTU 


PSAT/TPO 


TSAT/EBOIL DTSH/EBPL 


POUT /THL , 
RTDP/BBOIL ) 


C BEGIN INITIALIZATION FOR WATER BOILER 
C FLUID PROPERTIES - WATER BOILER DESIGN VALUES 
CPP = 1.0 
CPV = 0.474 
CH = 1.0 
HV = 875. 

USH = 108. 

UPH = 630. 

UBL = 880. 

WPDG = 45. 

WHDG = 770. 

AGAS = 8520.0 
BGAS = 15.426 
RGAS =0.596 
TABS = 460. 

ROWRE = 0.190 
ROWL = 59.4 

C PRESSURE DROP PARAMETERS 
AVPRS = 9.7 
R PH = 0.78 
R PS = 150. 

RBSP = 1.0 

RPHOP = ROWRE-RPS/ROWL 
C ITERATION PARAMETERS 
KCMAX = 500 
DK = 1.0 
ERMIN = 0.001 

C BOUNDARY CONDITIONS AT INLET - WP SPECIFIED 
C BOUNDARY CONDITION AT OUTLET - CHOKED NOZZLE 


BLCNOOOO 
BLCN0010 
BLCN0020 
BLCN0030 
BLCN0040 
BLCN0050 
B LCN0060 
BLCN0070 
BLCN0080 
BLCN0090 
BLCN0100 
BLCN0110 
BLCN0120 
BLCN0130 
BLCN0140 
BLCN0150 
BLCN0160 
BLCN0170 
BLCN0180 
BLCN0190 
BLCN0200 
BLCN0210 
BLCN022 0 
BLCN0230 
BLCN0240 
BLCN0250 
) BLCN0260 
BLCN0270 
BLCN0280 
BLCN0290 
BLCN0300 
BLCN03 10 
BLCN0320 
BLCN0330 
BLCN0340 
BLCN0350 
BLCN0360 
BLCN0370 
BLCN0380 
BLCN0390 
BLCN0400 
BLCN0410 
BLCN0420 
BLCN0430 
BLCN0440 
BLCN0450 
BLCN0460 
BLCN047 0 
BLCN0480 
B LCN0490 
BLCN0500 
BLCN0510 
BLCN0520 
BLCN0530 
BLCN0540 
BLCN0550 
BLCN0560 
BLCN0570 
BLCN0580 
BLCN0590 
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ICASE = 0 
TNK = 0.455 
C BOILER GEOMETRY 
ALT = 10. 

EPLUG = 0.215 
DBL = 0.049 
DPH = 0.049 
TUBNU =1.0 

C END INITIALIZATION FOR WATER BOILER 
C BEGIN GENERALIZED INITIALIZATIONS 

C INTERMEDIATE CONSTANTS - NTU VALUES AND RATIOS OF HEAT CAPACITY FLOWS 
A 1 = WPDG*CPV/(CH*WHDG) 

A2 = <HV*WPDG)/<CH*WHDG) 

A3 = (CPP*WPDG)/ (WHDG*CH) 

FNBL = 3.1416#TUBNU*UBL*DBL*ALT/(CH*WHDG) 

FNPH =3.1416*TUBNU*UPH*DPH*ALT/ (CPP*WPDG> 

FNPL = FNBL 

FNSH = FNBL*USH/ ( A1*UBL ) 

C END GENERALIZED INITIALIZATIONS 

C PREPARE FOR ITERATIONS AND READ INPUT VARIABLES 
1 CONTINUE 

READ (5*30 ) ICASE,KCMAX,DK,KPR 

WR I TE ( 6 1 33 ) ICASE, KCMAX , DK , KPR 

READ (5,31) WPO, PSATO, WH,TPO,THL,PO 

WRITE! 6,32) WPO, PSATO, WH,TPO,THL>PO 

WP =WPO/WPDG 

WH = WH/ WHDG 

PS AT = PSAT 0 

ASH =0.0 

DTSH =0.0 

PIN = PSATO 

WRITE (6,35) 

WRITE (6,34) WH,WP,RPH,RPS,UPH,USH,UBL, ALT, EPLUG 
WRITE (6,34) FNBL , FNPH, FNSH, R BSP, TNK, AVPRS,A1 ,A2, A3 
C MAIN BODY OF PROGRAM - VARY WORKING FLUID FLOW 
DO 54 KRW=1 ,25 
WP = AWP(KRW) 

C INITIALIZE FOR ITERATION 
K P RC = KPR 
KCNTR =0 
QUAL = 1.0 
DTSH = 0.0 
PSAT = PSATO * WP 
PIN = PSAT 

C BEGINNING OF ITERATION LOOP TO ADJUST OUTLET PRESSURE TO CORRESPOND 
C TO SPECIFIED FLOW. 

10 CONTINUE 

C COMPUTE HEAT TRANSFER PARAMETERS 
ASH = WP*A1/WH 
BNTU = FNBL/ WH 
PHNTU = FNPH/ WP 
APH = A3*WP/WH 
122 CONTINUE 

TSAT=AGAS/(BGAS-ALOG(PSAT) ) - TABS 
THLSH = THL - ASH*DT$H 

BBOIL =(A2*WP/(WH*(THLSH-TSAT) ) )*QUAL 
C TEST WHETHER AVAILABLE POWER IS ENOUGH TO VAPORIZE LIQUID. 

1 F ( BBO I L—O .99 >11,12,12 
12 CONTINUE 

C HEATING ENERGY LOW, REDUCE QUALITY AT OUTLET 


BLCN0600 

BLCN06 10 

BLCN0620 

BLCN0630 

BLCN0640 

BLCN0650 

BLCN0660 

BLCN0670 

BLCN0680 

BLCN0690 

BLCN0700 

BLCN07 10 

BLCN0720 

BLCN0730 

BLCN0740 

BLCN0750 

BLCN0760 

BLCN0770 

BLCN0780 

BLCN0790 

BLCN0800 

BLCN0810 

BLCN0820 

BLCN0830 

BLCN0840 

BLCN0850 

BLCN0860 

BLCN0870 

BLCN0880 

BLCN0890 

BLCN0900 

BLCN09 10 

BLCN0920 

BLCN0930 

BLCN0940 

BLCN0950 

BLCN0960 

BLCN097 0 

BLCN0980 

BLCN0990 

BLCN1000 

BLCN1010 

BLCN1020 

BLCN1030 

BLCN1040 

BLCN1050 

BLCN1060 

BLCN1070 

BLCN1080 

BLCN1090 

BLCN1100 

BLCN1110 

BLCN1120 

BLCN1130 

BLCN1140 

BLCN1150 

BLCN1160 

BLCN1170 

BLCN1180 

BLCN1190 
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o o 


DTSH = 0.0 
QUAL = QUAL - 0.001 
PSAT = 0 • 99 9* PS AT 
GO TO 122 

HEATING ENERGY SUFFICIENT 

COMPUTE LENGTHS OF REGIONS AND NTU OF PLUG 
11 CONTINUE 

BPHT =(TSAT-TP0)/(THLSH-(A2*WP/WH)*QUAL - TPO ) 

EL PH = AL0G(ABS(<1.0-APH)/(1.0-BPHT) ) )/<PHNTU*( l.O-APH)) 

IF (QUAL —1 .0 ) 112*113 *112 

112 CONTINUE 

ELSH = -0.00001 
EBOIL = 1.0 - ELPH 
GO TO 111 

113 CONTINUE 

RLBL = A LOG ( 1 .0/ ( 1 *0— BBOI L ) ) 

EBOIL = RLBL/BNTU 
ELSH = l.O-ELPH-EBOIL 
111 CONTINUE 

EBPL =EPLUG -ELPH 

C PRINT VARIABLES EVERY KPR ITERATIONS 
I F ( KPRC ) 4 1 *41 *42 

41 KPRC=KPR 
KSPR= 0 
WRITE(6*106) 

WRITE (6*107) WP*QUAL * PSAT * POUT * DP * TPL * TS AT * RBLD * RTDP 
WR I TE ( 6* 108 ) P IN *WH*TPO*THL* ELPH* ELSH* EBOIL* EBPL* BBOI L 

42 KPRC=KPRC-1 

C TEST WHETHER PREHEAT REGION EXTENDS OUTSIDE PLUG 
IF(EBPL)13*13*14 

13 CONTINUE 

C PREHEAT REGION EXTENDS BEYOND PLUG - SET BOILING AND SUPERHEAT PLUG 
C PARAMETERS TO ZERO 
PLNTU = 0.0 
ESHP = 0.0 
GO TO 15 

14 CONTINUE 

ESHP = EBPL - EBOIL 

C TEST FOR SUPERHEAT LENGTH* ESHP IN PLUG 
I F ( ESHP > 114*114*115 
115 CONTINUE 

C BOILING REGION ENDS BEFORE END OF PLUG - CALCULATE PLUG NTU 
EBPL = EBOIL 
PLNTU = FNPL*EBPL/WH 
GO TO 15 

114 CONTINUE 

C NO SUPERHEAT IN PLUG - CALCUALATE REGULAR BOILING PLUG NTU. 

ESHP = 0.0 

PLNTU = FNPL*EBPL/WH 

15 CONTINUE 

C TEST FOR OPERATION WITH EXIT QUALITY LESS THAN ONE 

C ELSH LESS THAN ZERO -ADJUST EXIT QUALITY. ELSH = ZERO - EXIT QUALITY 
C HAS BEEN ADJUSTED* CONTINUE NORMAL CALCULATION. ELSH GREATER THAN 
C ZERO - SOME SUPERHEAT * CONTINUE NORMAL CALCULATION 
IF(ELSH)17*16*16 
17 CONTINUE 

C BOILER OPERATION WITH QUALITY LESS THAN ONE 
ELSH = 0.0 
ASH = 0.0 
KQUAL = 100 


BLCN1200 

BLCN1210 

BLCN1220 

BLCN1230 

BLCN1240 

BLCN1250 

BLCN1260 

BLCN1270 

BLCN1280 

BLCN1290 

BLCN1300 

BLCN13 10 

BLCN1320 

BLCN1330 

BLCN1340 

BLCN1350 

BLCN1360 

BLCN1370 

BLCN1380 

BLCN1390 

BLCN1400 

BLCN1410 

BLCN1420 

BLCN1430 

BLCN1440 

BLCN1450 

BLCN1460 

BLCN1470 

BLCN1480 

BLCN1490 

BLCN1500 

BLCN1510 

BLCN1520 

BLCN1530 

BLCN1540 

BLCN1550 

BLCN1560 

BLCN1570 

BLCN1580 

BLCN1590 

BLCN1600 

BLCN16 10 

BLCN1620 

BLCN1630 

BLCN1640 

BLCN1650 

BLCN1660 

BLCN1670 

BLCN1680 

BLCN1690 

BLCN1700 

BLCN1710 

BLCN1720 

BLCN1730 

BLCN1740 

BLCN1750 

BLCN1760 

BLCN1770 

BLCN1780 

BLCN1790 
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171 
C 
C 


28 


CONTINUE 

ITERATE ON OUAL FOR A MAXIMUM OF 100 TIMES UNTIL BOILING LENGTH 
MATCHES END OF TUBE 
EBOIL = 1.0 - ELPH 
OBBL = 1.0 - E X P ( 0 • 0— EBOIL-BNTU) 

TEST FOR CONVERGENCE 

IF < ABS(QBBL-BB0IL )-0.001 
CONTINUE 

QUAL = QUAL*QBBL/BBOIL 
KOUAL = KQUAL - 1 
IF (KQUAL) 22,22,172 


) 111,111,28 


1.0) 174,174, 173 


172 CONTINUE 

IF (QUAL ~ 

173 CONTINUE 
QUAL « 1.0 
GO TO 122 

174 CONTINUE 
THLSH = THL 

BBOIL =( A2*WP/ (WH*(THLSH-TSAT) ) )*QUAL 

BPHT =(TSAT-TP0)/(THLSH-(A2*WP/WH)*QUAL - TPO ) 

ELPH = ALOG(ABS< (1.0-APH)/(1.0-BPHT) ) )/<PHNTU*( l.O-APH)) 

GO TO 171 
16 CONTINUE 

C OPERATION WITH SUPERHEAT - NORMAL CALCULATION 
SHNTU = FNSH*ELSH /WP 

DTSH = ( THL-T SAT )*( 1.0- ( 1 .O-ASH ) *EXP ( SHNTU* ( ASH-1 *0) ) ) 

TPL = TSAT +DTSH 
GO TO 18 

C COMPUTE PRESSURE DROPS 
18 CONTINUE 

ROW =PSAT/(ROWRE*RGAS*(TSAT+TABS )) 

PVT = WP*WP/ROW 

RLBL = AL0G(1.0/(1. 0— BBO I L ) ) 

C COMPUTE RATIOS OF PRESSURE DROP FOR EACH BOILER SECTION 
RSHDP = 1.0-ELPH-EB0IL +E SHP* ( RPS- 1 . 0 ) 

P PHD = ELPH - EPLUG 

C TEST FOR PREHEAT LENGTH OUTSIDE PLUG 
IF ( PPHD ) 60,60,61 

60 CONTINUE 

C PREHEAT REGION INSIDE PLUG 
RPHDP = RPH*ELPH *ROW 
GO TO 62 

C PREHAT REGION OUTSIDE PLUG 

61 RPHDP = ROW* ( RPH*EPLUG+RPHOP*PPHD ) 

62 CONTINUE 

RB - ( 1 .O-BBOIL ) /BBO I L 

RPLDP =( RBSP*RB*(RPS-1.0 )*( EXP( P LNTU ) -P LNTU- 1 . 0 ) / BNTU )*GUAL 

RSPDP = (RBSP*( 1 .0-RB*RLBL )/BNTU)*QUAL 

RTDP = RPHDP + RPLDP +RSPDP+RSHDP 

DP = PVT*RTDP*AVPRS 

RD = RTD P*P VT 

RBLDP = RPLDP+RSPDP 

C COMPUTE INLET AND OUTLET PRESSURES BASED ON CALCULATED PRESSURE DROP 
PIN = PSAT +RPHDP*PVT*AVPRS 
POUT = PIN -DP 

C WRITE INTERMEDIATE RESULTS EVERY KPR ITERATIONS 
IF(KSPR) 43,43,44 
43 K S PR = 1 

102 FORMAT ( 50H PRESOUT 
1 50H B RRSPLG PRDR BSP 


PRESDROP 
PRDR BOIL 


PR VAP TOT 
PRDRSHT 


PDROP 


PHEAT 

> 


BLCN1800 
BLCN1810 
BLCN1820 
BLCN1830 
BLCN1840 
BLCN1850 
BLCN1860 
BLCN1870 
BLCN1880 
BLCN1890 
BLCN1900 
BLCN1910 
BLCN1920 
BLCN1930 
BLCN1940 
BLCN1950 
BLCN1960 
BLCN1970 
BLCN1980 
BLCN1990 
BLCN2000 
BLCN2010 
BLCN2020 
BLCN2030 
BLCN2040 
BLCN2050 
BLCN2060 
BLCN2070 
BLCN2080 
BLCN2090 
BLCN2100 
BLCN2110 
BLCN2120 
BLCN2130 
BLCN2140 
BLCN2150 
BLCN2160 
BLCN2170 
BLCN2180 
BLCN2190 
BLCN2200 
BLCN2210 
BLCN2220 
BLCN22 30 
BLCN2240 
BLCN2250 
BLCN2260 
BLCN2270 
BLCN2280 
BLCN2290 
BLCN2300 
BLCN2310 
BLCN2320 
BLCN2330 
BLCN2340 
BLCN2350 
BLCN2360 
BLCN2370 
BLCN2380 
BLCN2390 
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103 FORMAT ( 1H ,8( E12.5) ) 

WRITE (6,102) 

WRITE (6,103) POUT ,DP , PVT , RPHOP , RPLDP ,RSPDP,RBLDP,RSHDP 
44 CONTINUE 

C COMPUTE BOILER FLOW BASED ON CHOKED NOZZLE AND COMPARE TO 
C SPECIFIED FLOW WP 

WP PR = TNK*POUT/SQRT(TPL +TABS ) 

WPERR = WP*GUAL - WPPR 

C TEST FOR CONVERGENCE ERROR IN FLOW LESS THAN ERMIN 
IF (ABS(WPERR)-ERMIN ) 20,20,19 

19 CONTINUE 

C ERROR STILL UNACCEPTABLE - CONTINUE ITERATION IF KCMAX NOT EXCEEDED 
KCNTR = KCNTR +1 
IF(KCNTR-KCMAX) 21,21,22 

21 CONTINUE 

C CASE I FLOW WP IS SPECIFIED, ITERATE ON PSAT ,ICASE =0 
PSAT = PSAT+DK*WPERR 
PIN = PSAT +RPHDP*PVT*AVPRS 
GO TO 10 

22 CONTINUE 

C NO CONVERGENCE - KCMAX EXCEEDED OR KQUAL = 0 

104 FORMAT ( 52H NO CONVERGENCE -KCMAX EXCEEDED 
WRITE (6,104) 

GO TO 20 

C NORMAL OUTPUT FOR SUCCESSFUL CONVERGENCE 

20 CONTINUE 

105 FORMAT ( 7H0 AFTER ,I4,23H ITERATIONS, FOR CASE =,I4, 

150H THE FOLLOWING BOILER VARIABLES HAVE RESULTED ) 

WRITE(6,105) KCNTR, ICASE 
WRITE (6,37) 

WRITE (6,36) APH,ASH,PHNTU,BNTU,SHNTU,PLNTU,ROW 
WR I T E ( 6 , 106 ) 

WRITE(6,107) W P , QUAL , P S A T , POUT , DP , TP L , TSAT , DT SH , RT DP 
WRITE(6,108)PIN,WH,TPQ,THL,ELPH,ELSH,EB0IL,EBPL,BB0IL 
WRITE (6,102) 

WRITE (6,103) POUT, RD, P VT, RPHDP, RPLDP, RSP DP, RB LDP, RSHDP 
54 CONTINUE 

C GO TO READ NEXT SET OF DATA 
GO TO 1 
END 

$DAT A 


45. 


0 


70. 


500 10.0 

770. 


100 

140. 410. 


85. 


BLCN2400 
BLCN2410 
BLCN2420 
BLCN2430 
BLCN2440 
BLCN2450 
BLCN2460 
BLCN2470 
BLCN2480 
BLCN2490 
BLCN2500 
BLCN25 10 
BLCN2520 
BLCN2530 
BLCN2540 
BLCN2550 
BLCN2560 
BLCN2570 
BLCN2580 
BLCN2590 
BLCN2600 
) BLCN26 10 
BLCN2620 
BLCN2630 
BLCN2640 
BLCN2650 
BLCN2660 
BLCN2670 
BLCN2680 
BLCN2690 
BLCN2700 
BLCN2710 
BLCN2720 
BLCN2730 
BLCN2740 
BLCN2750 
BLCN2760 
BLCN2770 
BLCN2780 
BLCN2790 
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APPENDIX C 


DESCRIPTION OF WATER BOILER 

The boiler is single-pass, counterflow type with single tube-in-shell construction, as 
shown in figure 1(b). It uses pressurized liquid water as the heating fluid in the shell, 
and water as the boiling fluid inside the tube. The boiling tube is standard stainless - 
steel tubing 5/8 inch (1. 59 cm) in outside diameter with a 0. 035-inch (0. 089-cm) wall. 

The shell is standard stainless- steel tubing 7/8 inch (2. 22 cm) in outside diameter with 
a 0. 049-inch (0. 125-cm) wall. The length of the heated section of the boiler is 10 feet 
(3. 05 m). Inside the boiling tube is an insert, to induce swirl flow, which consists of two 
sections. At the boiler inlet is a plug 2. 15 feet (0. 655 m) long machined to provide a 
spiral annular passage is 0. 050 inch (0. 127 cm) high and 0. 360 inch (0. 915 cm) wide. 
From the end of the machined plug to the boiler exit, the tube contains a 0. 067-inch- 
(0. 17 -cm-) diameter wire helically coiled to approximately 3 turns per inch (1. 18 
turns/ cm) which is attached to the inner surface of the boiling tube. In this region, 
concentric with the axis of the boiler, there is also a standard stainless -steel tube 
0. 25 inch (0. 635 cm) in outside diameter to support thermocouples. At the boiler exit, 
a valve serves to simulate a choked-nozzle boundary condition. At the boiler entrance, 
another valve actuated by a feedback controller on inlet flow, provides for precise 
regulation of the working -fluid liquid flow. 

The inlet and outlet pressures of the working fluid were measured by two Bourdon 
tube gages, each with an estimated accuracy of ±0. 2 psi (±1380 N/m ) (error, < 1/2 per- 
cent). The heating- and working -fluid inlet and outlet temperatures were measured with 
Chromel-Alumel thermocouples with an expected accuracy of ±1. 5° F (±0. 83° K) (error, 

< 1 percent). The heating -fluid flow was measured by a turbine flowmeter with an ac- 
curacy of ±15 pounds per hour (±6. 8 kg/hr) (error, <2.5 percent). The working -fluid 
flow was also measured by a turbine flowmeter with an accuracy of ±2 pounds per hour 
(±0.91 kg/hr) (error, < 5 percent). 
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APPENDIX D 


DESCRIPTION OF MERCURY BOILER 

The SNAP-8 boiler is a single-pass, counterflow heat exchanger. The mercury 
working fluid is distributed to seven parallel tubes, heated by liquid-metal NaK (the 
eutectic of sodium and potassium) flowing in a shell. The boiler tubes are 30 feet (9. 14m) 
long and double-wall construction where the inside wall is of tantalum, and the outside 
wall is an elliptically shaped stainless -steel tube. The small space between the two tubes 
is filled with stagnant NaK. The double-wall construction provides material compatibility 
and higher reliability. The only effect this special construction had on the analysis was 
to increase the relative magnitude of the wall heat-transfer resistance. Spiral plug 
inserts extend 4. 5 feet (1. 37 m) inside each tube. Sixteen parallel spiral grooves are 
cut into the plug to improve heat transfer and dispersion of the liquid. Analytically, the 
only difference between this multipassage plug and the single -passage plug of the water 
boiler is that the flow resistance to the mercury is less (i. e. , the r g 4 parameter is 
36 as compared with 150 for the water boiler). Only this dimensionless parameter is 
significantly different in the two boilers. The other parameters at design conditions, as 
listed in table n, are similar. This similarity in dimensionless parameters indicates 
that the pressure-drop characteristics should also be similar. The pressure-drop 
parameters for the design point of the mercury boiler, listed in the table, were evaluated 
by the methods described in appendix B. Representative values of the NTU parameters 
were calculated by using average constant coefficients of heat transfer estimated from 
the design calcualtions for the boiler (not shown here). An average overall heat transfer 
of 500 Btu per hour per °F per square foot (2830 J/(sec)(°K)(m^) was assumed for both 
the preheat and boiling regions and was based on an effective heat-transfer tube diameter 
of 1 inch (2. 54 cm). The coefficient of heat transfer in the superheat region was based 
solely on the working -fluid coefficient for design vapor flow. This coefficient was 
multiplied by the factor h w ^ re /hp^ n = 1.5 (see ref. 12 for s/D ratio of 36 to 
account for the spiral wire insert). The resultant Up, neglecting the heating -fluid and 
wall resistance was 80 Btu per hour per °F per square foot (454 J/ (sec)(°K)(m^)). The 
rest of the parameter values listed in table n can be calculated directly from the input 
variables and from the constants enumerated in the initialization part of the computer 
program listed in appendix E. 
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SIBFTC 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


BLCP LIST 

THIS IS A FORTRAN IV PROGRAM FOR COMPUTING THE PRESSURE DROP 
OF A SINGLE-PASS COUNTERFLOW BOILER. 

THE PROGRAM IS GENERAL EXCEPT FOR THE INITIALIZATION SECTION AND 
THE DATA CARDS. THE NUMERICAL CONSTANTS ARE SPECIFIED IN ENGLISH 
UNITS , ALTHOUGH ANY CONSISTENT UNITS CAN BE EMPLOYED. 


BTU PSI DEG-F ) 
SEVEN TUBE MERCURY 
CONSTANT PRESSURE 
OF WORKING FLUID 
WHICH THE PRESSURE 


BOILER 


FLOW 

DROP 


( IN UNITS 
FUNCTION 


30 

31 

32 

33 

34 

35 

36 

37 


(ENGLISH UNITS ARE IN FT HR LB 
THE NUMERICAL VALUES ARE FOR THE 
THE EXIT BOUNDARY CONDITION IS A 
THE DATA - AWP - SPECIFIES VALUES 
OF FRACTION OF DESIGN FLOW ) FOR 
IS TO BE EVALUATED. 

DIMENSION AWP ( 25 ) 

DATA AWP/ 1.0, 0.02, .06,0.1 ,0. 2*0. 3,0. 4,0. 5,0. 6, 0.7, 0.8, 0.9, 1.1, 
11.2,1.3,1.4,1.6,1.8,2.0,0.95,1.05,1.15,1.25,1.35,1.5 / 

FORMAT (2112, 1F12.7, 112 ) 

FORMAT ( 6F 1 2 • 7 ) 

FORMAT ( 49H WPO PSATO WH 


PO 

PARAMETERS 


C 

C 

c 


125H THL 

FORMAT ( 20H1 INPUT 
1 6H DK= , E 1 2 . 5 , 

FORMAT ( 1 H , 9 ( E 1 2 ♦ 5 ) ) 

FORMAT ( 50H1 WH/FNBL PDI/FNPH 
1 60H UPH/TNK USH/AVPRS UBL/A1 

FORMAT (1H , 7 ( E 1 2 . 5 ) ) 

FORMAT ( 50H APH ASH 

1 60H SHNTU PLNTU ROW 

106 FORMAT ( 50H WP/PIN QUAL /WH 

1 60H DP/ELPH TPL/ELSH TSAT/EBOIL 

107 FORMAT (1H , 9(E12.5)) 

108 FORMAT (1H , 9CE12.5)) 

BEGIN INITIALIZATION FOR MERCURY BOILER 
INPUT DATA FLOWS GIVEN FOR SEVEN TUBES 
BOILER GEOMETRY - MERCURY 7 TUBE BOILER 


1H , 6F 1 2 . 5 ) 

,6HICASE= r , I5/7H 
6H K PR= ,15 


RPH/FNSH 

ALT/A2 


TP- 

KCMAX= ,15, 
) 

RPS/RBSP 

EPLUG/A3 


PHNTU 

PSAT/TPO 

RBLD/EBPL 


BNTU 


POUT/THL 
RTDP/BBO I L 


DESIGN VALUES 


B4_ C P 0 0 0 0 
BLCP0010 
BLCP0020 
BLCP0030 
BLCP0040 
BLCP0050 
BLCP0060 
BLCP0070 
BLCP0080 
BLCP0090 
BLCP0100 
BLCP0110 
BLCP0120 
BLCP0130 
BLCP0140 
BLCP0150 
BLCP0160 
BLCP0170 
BLCP0180 
BLCP0190 
BLCP0200 
BLCP02 1 0 
8LCP022 0 
BLCP0230 
BLCP0240 
BLCP0250 
) BLC P0260 
BLCP0270 
BLCP0280 
BLCP0290 
BLCP0300 
BLCP0310 
B LCP0320 
BLCP0330 


ALT = 30.0 

BLCP0340 

EPLUG = 0.15 

BLCP0350 

DPH = 0.083 

BL C PO 360 

DBL = 0.083 

BLCP0370 

TUBNU = 7.0 

BLCP0380 

CPP = .033 

BLCP0390 

CPV = .025 

BLCP0400 

CH = 0.21 

BLCP0410 

HV = 123.1 

BLCP0420 

USH = 80. 

BLCP0430 

UBL = 500. 

BLCP0440 

UPH = 500. 

BLCP0450 

WPDG = 11500.0 

BLCP0460 

WHDG = 48000.0 

BLCP0470 

AGAS = 12100.0 

BLCP0480 

BGAS = 13.5 

BLCP0490 

RGAS = 0.0535 

BLC PO 500 

TABS = 460. 

BLCP0510 

ROWRE = 3.6 

BLCP0520 

ROWL = 760. 

BLCP0530 

= SSURE DROP PARAMETERS 

BLC PO 540 

AVPRS = 32.4 

BLCP0550 

R PS = 36. 

BLCP0560 

RPH = 0.28 

BLCP0570 

RBSP = 1.0 

BLC PO 5 80 

RPHOP = ROWR£*RPS/ROWL 

BLCP0590 
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c 

ITERATION PARAMETERS 

BLCP0600 


KCMAX = 500 

BLCP0610 


DK = 1.0 

BLCP0620 


ERMIN = 0.001 

BLCP0630 


POUTD = 270. 

BLCP0640 


PSMIN = 1.0 

BLCP0650 

c 

BOUNDARY CONDITIONS AT INLET - WP SPECIFIED 

BLCP0660 

c 

BOUNDARY CONDITION AT BOILER EXIT - CONSTANT PRESSURE 

BLCP0670 


ICASE = 0 

BLCP0680 

c 

END INITIALIZATION FOR MERCURY BOILER 

BLCP0690 

c 

BEGIN GENERALIZED INITIALIZATIONS 

BLCP0700 

c 

INTERMEDIATE CONSTANTS - NTU VALUES AND RATIOS OF HEAT CAPACITY 

FLOWS BLCP0710 


A 1 = WPDG*CPV/(CH*WHDG> 

BLCP0720 


A2 = (HV*WPDG)/(CH*WHDG) 

BLCP0730 


A3 = (CPP*WPDG)/(WHDG*CH) 

BLCP0740 


FNPH = 3.1416 5jc TUBNU*UPH J ! c DPH*ALT/ ( CPP*WPDG ) 

BLCP0750 


FNBL = 3.1416*TUBNU*UBL*DBL*ALT/(CH*WHDG> 

BLCP0760 


FNPL = FNBL 

BLCP0770 


FNSH = FNBL*USH/ ( A1*UBL ) 

BLCP0780 

c 

END GENERALIZED INITIALIZATIONS 

BLCP0790 

c 

PREPARE FOR ITERATIONS AND READ INPUT VARIABLES 

BLCP0800 

1 

CONTINUE 

BLCP0810 


READ ( 5 ? 30 ) I CASE, KCMAX ?DK?KPR 

BLCP0820 


WRITE(6?33) ICASE? KCMAX ? DK ? KPR 

BLCP0830 


READ ( 5 ? 3 1 ) WPO,PSATO?WH?TPO?THL?PO 

BLCP0840 


WRITE* 6? 32) WPO?PSATO?WH?TPO?THL?PO 

BLCP0850 


WP =WPO/WPDG 

BLCP0860 


WH = WH/ WHDG 

BLCP0870 


PSAT=PSATO 

BLCP0880 


ASH =0.0 

BLCP0890 


DTSH =0.0 

BLCP0900 


PIN = PSATO 

BLCP0910 


PD I = 0.0 

BLCP0920 


WRITE (6*35) 

BLCP0930 


WRITE ( 6? 34) WH? PD I ? RPH , R PS ? UP H ? USH , UBL ? AL T ? EP LUG 

BLCP0940 


WRITE ( 6? 34) FNBL ? FNPH? FNSH? RBSP ? TNK ? AVPRS? A1 ? A2? A3 

BLCP0950 

c 

MAIN BODY OF PROGRAM - VARY WORKING FLUID FLOW 

BLCP0960 


DO 54 KR W= 1 ? 2 5 

BLCP0970 


WP = AWP(KRW) 

BLCP0980 

c 

INITIALIZE FOR ITERATION 

BLCP0990 


K PRC = KPR 

BLCP1000 


KCNTR =0 

BLCP1010 


OUAL = 1.0 

BLCP1020 


DTSH = 0.0 

B LC P 1030 


PSAT = POUTD +60. 

B LC P 1040 


PIN = PSAT 

BLCP1050 

c 

BEGINNING OF ITERATION LOOP TO ADJUST SATURATION PRESSURE TO 

BLCP1060 

c 

CORRESPOND TO SPECIFIED FLOW 

BLCP1070 


10 CONTINUE 

B LC P 1 080 

c 

COMPUTE HEAT TRANSFER 

BLCP1090 


ASH = WP*A1/WH 

BLCP1100 


BNTU = FNBL/ WH 

BLCP1110 


PHNTU = FNPH/ WP 

BLCP1120 


APH = A3*WP/WH 

BLCP1130 

122 CONTINUE 

B LC P 1 1 40 


TSAT=AGAS/(BGAS-ALOG(PSAT) ) - TABS 

B LC P 1 1 50 


THLSH = THL - ASH*DTSH 

BLCP1160 


BBOIl =( A2*WP/<WH*(THl_SH-TSAT) ) )*OUAL 

BLCP1170 

C 

TEST WHETHER AVAILABLE POWER IS ENOUGH TO VAPORIZE LIQUID. 

BLCP1 180 


IF(BB0IL-0. 99)11, I2»12 

BLCP1190 
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o o 


12 CONTINUE 

C HEATING ENERGY LOW* REDUCE QUALITY AT OUTLET 
DTSH = 0*0 
QUAL = DUAL - 0.001 
GO TO 122 

HEATING ENERGY SUFFICIENT 

COMPUTE LENGTHS OF REGIONS AND NTU OF PLUG 
11 CONTINUE 

BPHT =(TSAT-TPO)/ < THLSH-I A2*WP/WH )*QUAL - TPO > 

EL PH = A LOG ( ABS ( ( 1 • O-APH ) / ( 1 . O-BPHT ) ) ) / ( PHNTU* ( 1. O-APH)) 

IF (QUAL -1.0) 112*113*12 

112 CONTINUE 

ELSH = -0.00001 
EBOIL = 1.0 - ELPH 
GO TO 111 

113 CONTINUE 

RLBL = ALOG< 1.0/ ( l.O-BBOIL ) ) 

EBOIL = RLBL/BNTU 
ELSH = l.O-ELPH-EBOIL 
111 CONTINUE 

EBPL =E PLUG -ELPH 

C PRINT VARIABLES EVERY KPR ITERATIONS 
I F ( KPRC ) 41 *41 *42 
41 KPRC = K PR 

K S PR= 0 


WRITE(6*106) 

WRITE ( 6*107) WP * QUAL * P SAT * POUT * DP * TPL * TS AT * RBLD * RTDP 
WRITE( 6,108 ) PIN, WH,TPO*THL* ELPH* ELSH* EBOIL, EBPL* BBO I L 
42 KPRC=KPRC— 1 

C TEST WHETHER PREHEAT REGION EXTENDS OUTSIDE PLUG 
I F ( EBPL ) 13*13*14 

13 CONTINUE 

C PREHEAT REGION EXTENDS BEYOND PLUG - SET BOILING AND SUPERHEAT PLUG 
C PARAMETERS TO ZERO 
ESHP = 0.0 
PLNTU = 0.0 
GO TO 15 

14 CONTINUE 

ESHP =EBPL - EBOIL 

C TEST FOR SUPERHEAT LENGTH* ESHP IN PLUG 
I F ( ESHP ) 114*114*115 

115 CONTINUE 

C BOILING REGION ENDS BEFORE END OF PLUG - CALCULATE PLUG NTU 
EBPL = EBOIL 
PLNTU = FNPL*EBPL/WH 
GO TO 15 
114 CONTINUE 

C NO SUPERHEAT IN PLUG - CALCUALATE REGULAR BOILING PLUG NTU. 

ESHP = 0.0 

PLNTU = FNPL*EBPL/WH 

15 CONTINUE 

C TEST FOR OPERATION WITH EXIT QUALITY LESS THAN ONE 

C ELSH LESS THAN ZERO -ADJUST EXIT QUALITY. ELSH = ZERO - EXIT QUALITY 
C HAS BEEN ADJUSTED* CONTINUE NORMAL CALCULATION. ELSH GREATER THAN 
C ZERO - SOME SUPERHEAT * CONTINUE NORMAL CALCULATION 
IF(ELSH)17, 16*16 
17 CONTINUE 

C BOILER OPERATION WITH QUALITY LESS THAN ONE 


ELSH 
ASH = 


« 0.0 
0.0 


BLCP1200 

BLCP1210 

BLCP1220 

BLCP1230 

BLCP1240 

BLCP1250 

BLCP1260 

BLCP1270 

BLCP1280 

BLCP1290 

BLCP1300 

BLCP1310 

BLCP1320 

BLCP1330 

BLCP1340 

BLCP1350 

BLCP1360 

BLCP1370 

BLCP1380 

BLCP1390 

BLCP1400 

BLCP1410 

BLCP1420 

BLCP1430 

BLCP1440 

BLCP1450 

BLCP1460 

BLCP1470 

BLCP1480 

BLCP1490 

BLCP1500 

BLCP1510 

BLCP1520 

BLCP1530 

BLCP1540 

BLCP1550 

BLCP1560 

BLCP1570 

BLCP1580 

BLCP1590 

BLCP1600 

BLCP1610 

BLCP1620 

BLCP1630 

BLCP1640 

BLCP1650 

BLCP1660 

BLCP1670 

BLCP1680 

BLCP1690 

BLCP1700 

BLCP1710 

BLCP1720 

BLCP1730 

BLCP1740 

BLCP1750 

BLCP1760 

BLCP1770 

BLCP1780 

BLCP1790 
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KQUAL = 100 

C ITERATE ON OUAL FOR A MAXIMUM OF 100 TIMES UNTIL BOILING LENGTH 
C MATCHES END OF TUBE 

171 CONTINUE 

EBOIL = 1.0 - ELPH 

OBBL = 1.0 - EXP ( 0 • 0- EBOIL-BNTU) 

C TEST FOR CONVERGENCE 

IF { ABS(QBBL-BBOIL J-0.001 ) 111*111*28 

28 CONTINUE 

QUAL = OUAL*OBBL/BBO I L 
KQUAL = KQUAL - 1 
IF (KQUAL) 22*22*172 

172 CONTINUE 

IF (QUAL - 1.0) 174*174*173 

173 CONTINUE 
QUAL = 1.0 
GO TO 122 

174 CONTINUE 
THLSH = THL 

BBOIL =(A2*WP/(WH*(THLSH-TSAT) ) )*QUAL 

BPHT ={TSAT-TP0)/(THLSH-(A2*WP/WH)*QUAL ~ TPO ) 

ELPH = ALOG(ABS< (1.0-APH)/(1.0-BPHT) )>/<PHNTU*( 1.0-APH)) 

GO TO 171 
16 CONTINUE 

C OPERATION WITH SUPERHEAT - NORMAL CALCULATION 
SHNTU = FNSH*ELSH /WP 

DTSH = (THL-TSAT >*( 1.0-( 1. O-ASH )*EXP(SHNTU*(ASH-1.0) ) ) 

TPL = TSAT + DTSH 
GO TO 18 

C COMPUTE PRESSURE DROPS 
18 CONTINUE 

ROW =PSAT/ ( ROWRE*RGAS*(TSAT+TABS )) 

PVT = WP*WP/ROW 

RLBL = ALOG( 1 .0/ ( 1.0— BBOIL ) ) 

C COMPUTE RATIOS OF PRESSURE DROP FOR EACH BOILER SECTION 
RSHDP = l.O-ELPH-EBOIL +E SHP* < RP S- 1 . 0 ) 

P PHD = ELPH - EPLUG 

C TEST FOR PREHEAT LENGTH OUTSIDE PLUG 
IF ( PPHD ) 60*60*61 

60 CONTINUE 

C PREHEAT REGION INSIDE PLUG 

RPHDP = (RPH-ELPH + PDI) * ROW 
GO TO 62 

C PREHAT REGION OUTSIDE PLUG 

61 RPHDP = ROW*( RPH*EPLUG +RPHOP*PPHD +PDI) 

62 CONTINUE 

RB =( 1 .O-BBOIL )/BBOIL 

RPLDP =(RBSP*RB*(RPS-1.0)*(EXP( P LNTU ) -P LNT U- 1 . 0 ) / B NTU )*QUAL 

RSPDP =( RBSP*( 1 .0-RB*RLBL ) /BNTU )*QUAL 

RTDP = RPHDP + RPLDP +RSPDP+RSHDP 

DP = PVT^RTDP^AVPRS 

RD = RTDP* PVT 

RBLDP = RPLDP+RSPDP 

C COMPUTE PIN AND PSAT BASED ON CALCUALTED PRESSURE DROP 
PIN = POUT D + DP 
PSATR = PSAT 

PSAT = PIN - RPHDP*PVT*AVPRS 
PER = PSAT - PSATR 
POUT = POUTD 

102 FORMAT ( 50H PRESOUT PRESDROP PR VAP TOT PDROP PHEAT 


BLCP1800 

BLCP1810 

BLCP1820 

BLCP1830 

BLCP1840 

BLCP1850 

BLCP1860 

BLCP1870 

BLCP1880 

BLCP1890 

BLCP1900 

BLCP1910 

BLCP1920 

BLCP1930 

BLCP1940 

BLCP1950 

BLCP1960 

BLCP1970 

BLCP1980 

BLCP1990 

BLCP2000 

BLCP2010 

BLCP2020 

BLCP2030 

BLCP2040 

BLCP2050 

BLCP2060 

BLCP2070 

BLCP2080 

BLCP2090 

BLCP2100 

BLCP2110 

BLCP2120 

BLCP2130 

BLCP2140 

BLCP2150 

BLCP2160 

BLCP2170 

BLCP2180 

BLCP2190 

BLCP2200 

BLCP22 1 0 

BLCP2220 

BLCP2230 

BLCP2240 

BLCP2250 

BLCP2260 

BLCP2270 

BLCP2280 

BLCP2290 

BLCP2300 

BLCP2310 

BLCP2320 

B LC P 2 330 

BLCP2340 

BLCP2350 

BLCP2360 

BLCP2370 

BLCP2380 

BLCP2390 
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1 50H B RRSPLG PRDR BSP PRDRBOIL PRDRSHT > 

C WRITE INTERMEDIATE RESULTS EVERY KPR ITERATIONS 
IF(KSPR) 43*43,44 

43 KS PR = I 
WRITE (6*102) 

103 FORMAT ( 1H *8(E12.5)) 

WRITE ( 6*103 ) POUT*DP*PVT*RPHDP*RPLDP*RSPDP*RBLDP* RSHDP 

44 CONTINUE 

C TEST FOR CONVERGENCE - ERROR IN PSAT LESS THAN 0.5 PSI 
IF < ABS < PER ) - 0.5 ) 20*20*45 

45 CONTINUE 

C ERROR STILL UNACCEPTABLE - CONTINUE ITERATION IF KC MAX NOT EXCEEDED 
PSAT = PSATR + PER/2.0 
KCNTR = KCNTR +1 
IF(KCNTR -KCMAX) 21*21*22 

21 CONTINUE 
GO TO 10 

22 CONTINUE 

C NO CONVERGENCE - KCMAX EXCEEDED OR KQUAL = 0 

104 FORMAT ( 52H NO CONVERGENCE -KCMAX EXCEEDED 
WRITE (6*104) 

GO TO 20 
20 CONTINUE 

C NORMAL OUTPUT FOR SUCCESSFUL CONVERGENCE 

105 FORMAT (7H0 AFTER *I4*23H ITERATIONS* FOR CASE =*I4* 

1 50H THE FOLLOWING BOILER VARIABLES HAVE RESULTED ) 

WR I TE ( 6 * 105 ) KCNTR* ICASE 
WRITE (6,37) 

WRITE (6*36) APH*ASH*PHNTU*BNTU*SHNTU*PLNTU*ROW 
WR I TE ( 6* 106 ) 

WRITE (6*107) WP *OUAL * P S AT * POUT * DP * TPL * TS AT * DT SH * RTDP 
WRITE (6* 1 08 ) P IN* WH*TPO*THL* EL PH* EL SH,EBOIL* EBPL*BBOIL 
WRITE (6*102) 

WRITE (6*103) POUT * RD * P VT * RP HDP * RP LDP *RSPDP*RBL DP* RSHDP 
54 CONTINUE 

C GO TO READ NEXT SET OF DATA 



GO TO 1 
END 





$DATA 

0 

500 10.0 

100 



11500 

. 

270. 48000. 

500. 

1300. 

360. 


BLCP2400 
BLCP2410 
BLCP2420 
BLCP2430 
BLCP2440 
BLCP2450 
BLCP2460 
BLCP2470 
BLCP2480 
BLCP2490 
BLCP2500 
BLCP2510 
BLCP2520 
BLCP2530 
BLCP2540 
BLCP2550 
BLCP2560 
BLCP2570 
BLCP2580 
) BLC P2 590 
BLCP2600 
BLCP2610 
BLCP2620 
BLCP2630 
BLCP2640 
BLCP2650 
BLCP2660 
BLCP2670 
BLCP2680 
BLCP2690 
BLCP2700 
BLCP2710 
BLCP2720 
BLCP2730 
BLCP2740 
BLCP2750 
BLCP2760 
BLCP2770 
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TABLE I. - PARAMETERS AND INPUT VARIABLES FOR DESIGN POINT OF WATER BOILER 


(a) Parameters 


Dimensionless parameter 

Definition 

Design 

value 

Effectiveness of preheat region, €p 

( T sat - T wi)/( T hb - T wi) 

0. 823 

Effectiveness of boiling region, 

X 0 H w w w / ( T hs - T sat) w h c h 

0. 575 

Number of transfer units (NTU) of boiler based on overall 
average heat -transfer coefficient in preheat region, Np 

VVt /w w c w,Z 

21.6 

NTU of boiler based on overall average heat -transfer 

tWt/'Vh 

1.76 

coefficient in boiling region, 



NTU of boiler based on overall average heat-transfer 
coefficient in superheat region, N g 

U s’ rD sV w w c w,v 

3. 36 

Ratio of working -fluid to heating -fluid heat capacity in 

W w C w,Z /w h c h 

0. 0584 

preheat region, ap 


Ratio of working -fluid to heating -fluid heat capacity in 


0. 0277 

super -heated region, a g 


Preheat length (dimensionless), Lp 

Vh 

0. 083 

Superheat length (dimensionless, L g 

Vh 

0. 431 

Boiling length (dimensionless), L^ 

l ^ l t 

0. 485 

Plug length, L pl 

z pi^t 

0. 215 

Ratio of preheat -plug -region pressure gradient to super- 

V f s )( fv ) d/%,p(p p A)( A s/V (dz pi /dz) 

0.78 

heat pressure gradient, r 

p, b 


Ratio of boiling to superheat friction factors, r^ g 

^b, s^s 

1. 0 

Ratio of pressure gradients for vapor at rated flow in plug 

( f P i,b/ f b,sHV p s )( VV (dz pi /dz) 

150. 0 

region to those of vapor at rated flow in swirl -wire 
region, r g ^ ^ 



(b) Input variables 


Variables 

Design value 

Normalizing pressure-drop factor based on pressure gradient in superheat region at 

9.7 (67 000) 

design flow, AP y psia (N/m^ abs) 


Working -fluid outlet pressure, P Q psia (N/m^ abs) 

64. 5 (445 0C0) 

Working -fluid design flow, w w lb/hr (kg/hr) 

45 (20. 4) 

Heating -fluid design flow, w^ lb/hr (kg/hr) 

770 (349. 0) 

Working -fluid inlet temperature, T wi °F (°K) 

140 (333) 

Heating -fluid inlet temperature, T^ d , °F (°K) 

410 (483) 

Total boiler tube length, Z^., ft (m) 

10 (3. 048) 

3 3 

Working-fluid vapor density at design flow in superheat region, p ,, lb/ft (kg/m ) 

V y U 

0. 155 (2. 48) 

Working -fluid liquid density, p 7 , lb/ft ^ (kg/m^) 

w, t 

59. 4 (950) 

Hydraulic diameter D in the boiling, preheat, and superheat regions, ft (m) 

0. 049 (0. 015) 

Working -fluid saturation temperature, T °F (°K) 

320 (433) 

f 9 

Working -fluid saturation pressure, P ga ^ ^ , psia (N/m ) 

89 (615 000) 
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TABLE E. - COMPARISON OF MERCURY AND WATER BOILERS 


PARAMETERS AND DESIGN INPUT VARIABLES 
(a) Parameters 


Parameter 

Mercury 

Water boiler 


boiler 

(1 of 7 tubes) 



Desigi 

value 

Effectiveness of preheat region, 

0.908 

0. 823 

Effectiveness of boiling region, 

.698 

0. 575 

Number of transfer units (NTU) of boiler based on overall average 

72. 1 

21.6 

heat-transfer coefficient in preheat region, N p 



NTU of boiler based on overall average heat-transfer coefficient 

2.72 

1.76 

in boiling region, N fo 



NTU of boiler based on overall average heat-transfer coefficient 

15.2 

3.36 

in superheat region, N g 



Ratio of working -fluid to heating -fluid heat capacity in preheat 

0. 0376 

0. 0584 

region, a p 



Ratio of working -fluid to heating -fluid heat capacity in super- 

0. 0285 

0. 0277 

heated region, a g 



Preheat length (dimensionless), L p 

0. 0336 

0. 083 

Superheat length (dimensionless), L g 

0. 525 ! 

0. 431 

Boiling length (dimensionless), L^ 

0. 441 

0. 485 

Plug length, L pl 

0. 15 

0.215 

Ratio of preheat -plug-region pressure gradient to superheat 

0.28 

0.78 

pressure gradient, r, e 

p, b 



Ratio of pressure gradients for vapor at rated flow in plug region 

36.0 

150.0 

to those of vapor at rated flow in swirl -wire region, r g ^ 



Ratio of boiling to superheat friction factors, r fe g 

1.0 

1.0 


(b) Input variables 


Input variable 

Mercury boiler 
(1 of 7 tubes) 

Water boiler 


Design value 

Working -fluid inlet temperature, T^ d , °F (°K) 

500 (533) 

140 (333) 

Heating -fluid inlet temperature, T hi d , °F (°K) 

1300 (978) 

410 (483) 

Working -fluid flow, w w d , lb/hr (kg/hr) 

1644 (747) 

45 (20. 4) 

Heating -fluid flow d , lb/hr (kg/hr) 

6860 (3120) 

770 (349. 0) 

Working -fluid outlet pressure, P Q d , psia (N/m 2 abs) 

270 (186X10 6 ) 

64. 5 (445 000) 

Normalizing pressure-drop factor based on pressure gradient 

32.4 (2.24xl0 5 ) 

9.7 (67 000) 

in superheat region at design flow, AP y d , psia (N/m 2 abs) 



Total boiler tube length, Z^, ft (m) 

30 (9. 14) 

10 (3. 048) 

Working -fluid vapor density at design flow in superheat 

3.6 (57.7) 

0. 155 (2.48) 

region, p y d , lb/ft 3 (kg/m 3 ) 



Working -fluid liquid density, p w ^ lb/ft ^ (kg/m 3 ) 

760 (12 200) 

59. 4 (350) 

Working -fluid saturation temperature, °F (°K) 

1100 (866) 

319 (433) 

Hydraulic diameter D in preheat, boiling, and superheat 

0.083 (0.025) 

0. 049 (0. 015) 

regions, ft (m) 



_ 




Overall heat-transfer coefficient 


Temperature 




Figure 3. -Temperature and overall heat-transfer-coefficient profile 
showing difference in actual and assumed characteristic. 



Figure 4. - System representation of boiler equations for superheat at exit, 
at exit, choked nozzle. 


inlet flows, W w and w h ; specified inlet temperatures, T w j and T h j; boundary condition 


















Quality at exit Superheat at exit (T hi - T $at d = 90°F; Boiler outlet pressure, P 0 /P 0 d Pressure-drop function (o. 

3050 K ). (Two-W'ThrW ' ’ 






(d) Variation in quality at boiler exit. 


Figure 5 . - Comparison of theory with water-boiler data for variation in working-fluid flow. Working- 
fluid design flow, 45 pounds per hour (20. 4 kg/hr); ratio of pressure gradients for vapor at rated flow 
in plug region to vapor at rated flow in swirl -wire region, 150; ratio of preheat -plug-region pressure 
gradient to superheat pressure gradient, 0.78; ratio of boiling to superheat friction factors, 1.0; 
normalizing pressure-drop factor based on pressure gradient in superheat region at design flow, 

9. 7 psi (6. 7xl(r N/m 2 ); working-fluid inlet temperature (design value), 140° F (333° K); heat-fluid 
inlet temperature (design value), 410° F (483° K); heat -fluid design flow, 770 pounds per hour 
(349 kg/hr); choked nozzle at boiler exit set for 64 . 5 psia (445 000 N/m 2 abs) at design flow. 
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Pressure-drop function, <p t 


4 1 — 


O Experimental 
Calculated 



Variation of heating-fluid flow, W h , fraction of design 

(a) Variation in heating-fluid flow. Input variable: Heating fluid 
design flow, 45 pounds per hour (20.4 kg/hr); choked nozzle at exit. 
Theory parameters: ratio of preheat-plug-region pressure gradient 
to superheat pressure gradient, 0. 78; ratio of pressure gradients 
for vapor at rated flow in plug region to those of vapor at rated flow 
in swirl-wire region, 150; ratio of boiling to superheat friction 
factors, 1.0; normalizing pressure-drop factor based on pressure 
gradient in superheat region at design flow, 9. 7 psi (6. 7xl0 4 N/m 2 ). 



(b) Variation in heating-fluid inlet temperature. Working-fluid 
saturation temperature (design value), 320° F (433° K). 


3,— 



.4 .6 .8 1.0 1.2 1.4 1.6 


(c) Variation in working-fluid inlet temperature. Working-fluid 
saturation temperature (design value), 320° F (433° K). 

Figure 6. - Comparison of data with predicted pressure drop as 
function of input variables. Heating-fluid inlet temperature (de- 
sign value), 410° F (483° K); working-fluid inlet temperature (de- 
sign value), 140° F (333° K); dimensionless heating-fluid design 
flow, 1.0; dimensionless working-fluid design flow, 1.0. 
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Pressure-drop function, <p^ Dimensionless length, L Dimensionless pressure-drop function 


2.0 


Pressure-drop function 



(a) Tube with multipassage plug insert. Ratio of pressure gradients for vapor at rated flow in plug region 
to those of vapor at rated flow in swirl-wire region, 36; dimensionless plug length, 0. 15. 




(c) Tube without plug. Ratio of pressure gradients for vapor at rated flow in plug region to those of 
rated flow in swirl-wire region, 1.0. 


Figure 7. - Pressure drop and interface positions as function of working-fluid flow showing compo- 
nents due to superheat, boiling, and preheat regions. SNAP -8 mercury-boiler parameters: ratio 
of boiling to superheat friction factors, 1.0; ratio of preheat-plug-region pressure gradient to super- 
heat pressure gradient, 0. 28; normalizing pressure-drop factor based on pressure gradient in super- 
heat region at design flow, 32.4 psi (2.24xl0 5 M/m 2 ); heating-fluid design flow, 6860 pounds per hour 
(3120 kg/hr); working-fluid design flow, 1644 pounds per hour (747 kg/hr); working-fluid outlet pres- 
sure (design value), 270 psia (1. 86xl0 6 N/m? abs); heating-fluid inlet temperature (design value), 
1300° F (978° K) ; working-fluid inlet temperature (design value), 500° F (533° K); choked nozzle at 
exit; no orifice at inlet. 
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Dimensionless pressure drop, ^ + rj 



Figure 9. - Pressure drop as function of working-fluid flow showing effect of various inlet orifices 
{ Tj parameter) and difference between choked nozJe and constant exit pressure. Working-fluid 
design flow, 1649 pounds per hour (747 kg/hr); dimensionless heating-fluid design flow, 1.0; ratio 
of boiling to superheat friction factors, 1.0; ratio of preheat-plug-region pressure gradient to super- 
heat pressure gradient, 0. 28; ratio of pressure gradients for vapor at rated flow in plug region to 
those of vapor at rated flow in swirl -wire region, 36; normalizing pressure-drop factor based on 
pressure gradient in superheat region at design flow, 32.4 psia (2.24x10^ N/m^ abs); heating-fluid 
inlet temperature (design value), 1300° F (978° K) ; working-fluid inlet temperature (design value), 
500° F (533° K). 
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Dimensionless pressure drop, + r { 





Pressure-drop function, ^ Pressure-drop function, 


Dimensionless 
heating-fluid flow, 
w h /w M 


(a) Effect of heating-fluid flow variation. 





.2 .4 .6 .8 1.0 1.2 1.4 1.6 1 

Working-fluid flow, W w , fraction of design 

(b) Effect of heating-fluid inlet-temperature variation. Working-fluid saturation temperature (design 
value), 1100° F (866° K). 

Figure 11. - Boiler pressure-drop function as function of working-fluid flow showing effect of heating 
fluid parameters. Heating-fluid design flow, 6860 pounds per hour (3120 kg/hr); heating-fluid inlet 
temperature (design value), 1300° F (978° K) ; working-fluid inlet temperature (design value), 

500° F (533° K); working-fluid outlet pressure (design value), 270 psia (constant) (1.86x10° N/m?); 
working-fluid design flow, 1644 pounds per hour (747 kg/hr); constant exit pressure. 



Pressure-drop function 
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l cm) high by 0.360 in. (0.914 cm) 
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0.067-in. (0.17-cm) 
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inner wall of 3/8-in. 
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Heater-loop 
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Reference for 
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Working-loop 
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7/8-in.- (2.22-cm-) o.d. j / / 

tubing; wall thickness, j I / 
0.049 in. (0.1245 cm) _j j j 
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tubing; wall thickness, / / 

0.035 in. (0.089 cm) -/ j 
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tubing; wall thickness, 
0.049 in. (0.1245 cm) 




Tubethermo- 
= 3 — - couple leads 
outlet 

\ 

L 1/4-in. -(0.635-cm-) o.d. 
tubing; wall thickness, 
0.049 in. (0.1245 cm) 
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Figure 12. - Water-boiler cross section. 
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